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FOREWORD 


This report summarizes the major technical findings made in 
the research program at the University of Florida sponsored hy 
NASA Lewis Research Center under Grant No. NSG-3018. This grant 
is a continuation of the support began under the same grant number 
on Jvine 24, 1974. We report here our findings during the period 
of the grant: July 24, 1975 to July 24, 1976. 

The work to be reported came from the efforts of several key 
people. C. T. Sah suggested to NASA the possibility of placing this 
grant research at the University of Florida, and afterwards, together 
with the author, helped make the arrangements and formulate the tech- 
nical plan of the research. Dr. Sah has contributed significantly to 
nearly every aspect of the research described in this report. After 
joining the faculty of the University of Florida in September, 1976, 
Arnost Neugroschel became a key contributor to the research program. 

The technical collaboration of M. P. Godlewski and H. W. Brandhorst, 
Jr., of NASA Lewis Research Center greatly aided the completion of 
several phases of the research. During the tenure of the NASA grant, 
the Energy Research and Development Administration began supporting 
a separate but related research program at the University of Florida 
under the author's direction. This support aided the research reported 
here, and made possible a valuable and continuing technical interaction 
with J. G. Fossum of Sandia Laboratories (ERDA) , Albuquerque, New 
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Mexico. E. L. Burgess of Sandia tiaboratories also contributed to 
the research findings that we report. In addition to these senior 
Investigators, several graduate students at the University of Florida 
participated in the research. They are; P. J. Chen, S. C. Pao, 

G. Badeaux, M. P. Callaway, and E. Therez. 

Our research sponsored by NASA during the year (July 1975 to 
July 1976) produced the following technical papers: 

1. P. A. Lindholm and C. T. Sah, "Process- Induced Defects in Terrestrial 
Solar Cells," Proceedings of the First ERPA . Semiannual Solar Con- 
version Program Conference , Los Angeles, pp. 578-590, July 1975 
(invited paper) . 

2. M. P. Godlewski, H. W. Brandhorst, Jr., F, A. Lindholm and C. T. 

Sah, "Experimental Investigation of the Diffused Layer of a 0.1 
ohra-cm Silicon Solar Cell," Symposia of the 1976 Device Research 
Conference and the 1976 Electronic Materials Conference , Salt Lake 
City, Utah, Jme 23, 1976, . ' 

3. F. A. Lindholm and C. T, Sah, "Studies of Basic Mechanisms Influencing 
Solar Cell Efficiency for Terrestrial Applications," Proc. .National 
Solar Photovoltaic Program Review Meeting, session III, paper 7, 

Orono, ME, 1976. 

4. F. A. Lindholm, A. Neugroschel, and C. T. Sah, "Basic Mechanisms 
Governing Solar-Cells Efficiency," Proc. Int. Solar Energy Society , 
Vol. 6, pp, 120-129, Winnipeg, Manitoba, Canada, Aug. 1976. 

5. F, A. Lindholm and C. T. Sah, "Normal Modes of pn Junction Devices 
for Mate rial- Parameter Determination," J. Appl. Phys. , Vol. 47, pp. 
4203-4205, Sep. 1976. 

6. F- A. Lindholm, C. T. Sah,' A. Neugroschel, M. P. Godlewski, and 

H. W. Brandhorst, "Methodology for Experimental Determination of 
Gap Shrinkage and Lifetimes in Emitter and Base of pn Junction 
solar cells," Record of Twelfth Photovoltaic Specialists Conf . , 
paper 1.1, Nov. 1976. 

7. F. A. Lindholm, J. G. Fossum, and E. L. Burgess, "Basic Corrections 
to Solar Cell Performance Required by Inapplicability of Super- 
position," Record of Twelfth Photovoltaic Specialists Conf. , paper 

I. 6, Nov. 1976. 

8. F. A. Lindholm, A. Neugroschel, and C. T. Sah, "Junction Modeling 
for Solar Cells — Theory and E>^eriment," Tech. Digest 1976 Int. 
Electron Devices Meeting , 76CH1151-OED, pp. 61-65, Dec. 1976 
(invited paper) . 
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9. F. A. Lindholm, C. T. Sah, and J. G. FossuiU/ "Fundamental Considera- 
tions for Improving the Efficiency of Junction Silicon Solar Cells," 
Proc. National Solar Photovoltaic Program Review Meeting , session 
III, paper 1, San Diego, Calif., Jan. 1977. 

10. F. A. Lindholm and C, "i‘. Sah, "Fundamental Electronic Mechanisms 
Limiting the Performance of Solar Cells," IEEE Trans. Electron 
Devices , April 1977, 

11. C. T. Sah and F. A. Lindholm, "Carrier Generation, Recombination, 
Trapping and Transport in Semiconductors with Position-Dependent 
Composition," IEEE Trans. Electron Devices , April 1977. 

12. J. G. Fossum and F. A. Lindholm, "The Dependence of Open-Circuit 
Voltage on Illumination Level in p-n- Junction Solar Cells," IEEE 
Trans, on Electron Devices , April 1977, 

13. F. A. Lindholm, A. Neugroschel, C. T. Sah, M. P. Godlewski, and 

H. W. Brandhorst, Jr., "Methodology for the E:i^erimental Determina- 
tion of Gap Shrinkage and Lifetimes in the Emitter and Base of p-n 
Junction Solar Cells and Other p-n Junction Devices," IEEE Trans. 
Electron Devices , April 1977. 

14. A. Neugroschel, F. A. Lindholm, and C. T. Sah, "A Method for 
Determining the Emitter and Base Lifetimes in p-n Junction Diodes," 
IEEE Trans. Electron Devices, June 1977. 


15. M. P. Godlewski, H. W. Brandhorst, Jr., F. A. Lindholm, and C. T. 

Sah, "Ejjperimental Investigation of Excess Charge and Time Constant 
of Minority Carriers in the Diffused Layer of 0.1 ohm-cm Silicon 
Eolar Cell , " Int. J. Electronic Materials , accepted for publication . 

In addition, the author gave seminars on results of our NASA- sponsored 
research at: Oak Ridge National Laboratories, the University of North 

Carolina {Chapel Hill), Sandia Laboratories, Louisiana State University, 
Communication Satellite Corporation Laboratories, (Clarksburg, Maryland) and 
the Institute of Energy Conversion of the University of Delaware. 
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SUMMARY 


Single-crystal silicon p-n- junction solar cells made with low- 
resistivity substrates show poorer solar energy conv-er^iion efficiency 
than traditional theory predicts * The purpose of the research sponsored 
by grant NSG-3018 is to identify and characterize the physical mechanisms > 
responsible for this discrepancy. Attention concentrates on the open- 
circuit voltage in shallow- jiinction cells of about 0.1 ohm-cm sub- 
strate resistivity under AMO {one-stin) conditions. 

The open— circuit voltage depends on the electronic parameters that 
characterize the transport, recombination, and generation of electrons 
and holes in silicon. The many basic mechanisms that can control these 
electronic parameters are reviewed. Mechanisms contributing to high re- 
combination rates (low lifetime) and large minority-carrier storage 
(effective gap shrinkage) are identified as the two main origins of the 
discrepancy in a one-dimensional model of the solar cell. The in- 
adequacies of existing theoretical models and the lack of relevant ex- 
perimental data and of methods capable of yielding that data prevent 
an assessment of the relative importance of gap shrinkage and low life- 
time in any particular solar-cell structure. Such an assessment is 
needed to make possible the identification of the basic mechanisms con- 
trolling Thus .much of the effort reported here focusses on the 

development of experimentally-based procedures for determining the ef- 
fective lifetime and the minority-carrier storage in the quasi-neutral 
emitter. 

A method is described that provides an experimental means for the 
first time to separate and determine the emitter and base lifetimes in 
a p-n diode after the junction has been fabricated. In the method, 
several static and transient measurements are analyzed using physical 
models of the diode characteristics. To illustrate the method, diffused 
silicon diodes are fabricated having substrate (base) impurity concentra- 
tions ranging from 10 to nearly 10 phosphorous atoms per cubic centi- 
meter. The results show an emitter lifetime that is much smaller than 
the base lifetime in the diode having the highest base doping concentra- 
tion. In this diode, the recombination current from the emitter is 65% 
of the recombination current from the base, demonstrating the significance 
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of the emitter in governing the static current-voltage dependence. 

The importance of emitter recombination to the transient characteristics 
is also demonstrated- The paper emphasizes the techniques by which the 
base and emitter lifetimes are distinguished. It also demonstrates the 
need for carefully basing the quantitative analysis of the measiirements 
on the underlying diode physics. 

Based on the understanding gained from the study just described, 
an experimentally-based methodology is then developed that determines 
the effective gap shrinkage and lifetime in the emitter of a p-n- junc- 
tion solar cell. It provides the first experimental means available for 
assessing the importance of gap shrinkage relative to that of large re- 
combination rates in the highly-doped emitter. As an additional result 
of the procedures employed, the base lifetime is also determined. The 
methodology pertains to a solar cell after the junction is formed. 

Hence each material parameter determined includes the effects of the 
processing used in junction fabrication. The methodology consists of 
strategy and procedures for designing experiments and interpreting data 
consistently with the physical mechanisms governing device behavior- 
This careful linking to the device physics \incovers the material para- 
meters concealed in the data. To illustrate the procedures , they are 
applied to an n -p solar cell having a substrate resistivity of about 
0 , 1 ohm-cm. 

The principle of superposition is used to derive from fundamentals 
the widely used approximation that the current-voltage characteristic of 
an illuminated solar cell is the dark current- voltage characteristic 
shifted by the short-circuit current. Thus the derivation requires the 
linearity .of the boundary-value problems that underlie the electrical 
characteristics. This focus on linearity defines the conditions that 
must hold if the shifting approximation is to apply with good accvufacy. 

For the solar cell of interest (a 0-1 ohm-cm silicon cell under AMO 
conditions) , the shifting approximation is shown to be valid tinless material 
parameters such as the minority-aarrier lifetimes change appreciably under 
AMO sunlight from the values they have in the dark. Measurements of the 
dark and illuminated current-voltage characteristics of a 0.1 ohm-cm 
silicon cell indicate that this does not happen and directly validate 
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the shifting approximation for such a cell. The theoretical study des- 
cribed here suggests that the shifting approximation will likely be 
valid for all single-crystal silicon cells under low-injection condi- 
tions in the absence of sizable series resistance. 

Use of the shifting approximation enables the discrepancy in open- 
circuit voltage seen in single-crystal silicon cells under AMO conditions . 
to be investigated by studying the electrical behavior of the' cell in 
the dark. Thus it underlies many of the methods and findings- described 
in this report. 
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CHAPTER I, INTRODUCTION 

The solar energy conversion efficiency n of a solar cell is 
controlled by the product of three factors: the short-circuit photo- 

current (I_„) r the curve factor or fill factor (FF) , and the open- 
circuit phptovoltage Cv^^) . For high-efficiency silicon p-n-junctioh 
cells, the values of FF and I approach closely the limits- thought 
to be technologically achievable, but the values of occurring 

in commercially available silicon cells fall far short of theoretical expecta- 
tions [1] , Thus the greatest gains possible in n will likely result 

from increasing 

Figure 1 illustrates the discrepancy between theoretical and 
experimental V . As the substrate (or base) doping concentration 

Ng increases, the experimentally observed peaks at a substrate 

. . 17-3 

resxstivity of about 0.1 ohm-cm, corresponding to M 4 x 10 cm 

5 

for a p-type substrate. This maximiom V of about 600 mV falls 

OC 

nearly 100 mV below the value of V^^ predicted by classical theory 

for this base resistivity. Removal of this discrepancy in V^^ 

would thus increase n from the maximum of approximately 14% (AMO) 

17 —3 

now seen [1] to above 17%. Increases of N beyond N = 4 x 10 cm 

B B 

are unlikely to produce further improvements in n because the minor- 

18 “*3 

ity-carrier base lifetime falls sharply as N approaches 10 cm 

B 

The purpose of the research program supported by NASA Grant NSG-3018 
is to explain, in fundamental terms, the reasons for the discrepancy that 
exists between the observed values of V^^ and those predicted by traditional 
theory. The research program concentrates on V^^ in (n-on-p) single- 
crystal silicon p-n- junction solar cells. 
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1.1 Anomalous Values of Minority-Carrier Storage and of 
Effective Lifetimes in the Emitter _ 

To explore the origin of the discrepancy, note first that 
the widely-used shifting approximation , 

KV) 5 I^(V) - (1) 

provides for 0.1 ohm-cm silicon cells an accurate description under 
AMO conditions [2j . This approximation views the current-voltage 
characteristic I (V) of the illuminated cell as the dark characteristic 
Ij^(V) shifted by the short-circuit current The physical mech- 

anisms occurring in the dark that produce (V) thus become of cen- 
tral importance. 

By the principles of charge control, 


I^(V) 


IS. 






^B 


(2) 


Here Q„/t_ is the recombination current in the quasi-neutral emitter, 

£i £i 

described as the ratio of excess carrier charge to the phenomen- 
ological lifetime in the emitter, and is the corresponding 

recombination current in the quasi-neutral base. At voltages near 
. ^QQt ihe recombination current in thie junction space-charge region 
of silicon cells contributes negligibly to lg(V) [3-6]. 

The open-circuit voltage establishes the balance 


° - S'V' - "sc 


(3) 


Because I , Q , and Q all increase exponentially with increasing 
D B 

voltage V, the smaller is lg(V) then the larger will be Classical 

theory assumes that I = the emitter recombination current 

being assumed negligible because the relatively heavy doping concen- 
tration in the emitter is assumed to make Q for AMO conditions, however, 

£ B 
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measured values of V for a O.lfi-cm cell [1-4] imply that Q_/t » 

OC ' E S 

This anomalously large emitter recombination current occurs 

either because Q is larger than is predicted classically, or T_ 

E E 

is smaller, or for both of these reasons. Various physical mech- 
anisms have been suggested to justify the existence of an anoma- 
lously large minority carrier storage or an anomalously small 
lifetime in the emitter [7,8]. 

1. 2 Order of Presentation 

Chapter 2 sturanarizes the many fundamental electronic mechanisms 
that could contribute to the anomalous values of minority- carrier 
charge storage and lifetime. Chapter 3 concentrates on lifetime. 

It describes methods that enable for the first time the determina- 
tion by electrical measurements of the .phenomenological lifetime 

T in the emitter. The methods give also the value of the base 
E 

lifetime and identify and characterize the imperfection centers 

in the base- In Chapter 4, the methods are extended so that they 

yield not only the phenomenological lifetime but also the minority- 

carrier charge storage Q in the emitter. The application of these 

methods to an n'^-p junction solar cell having a base resistivity 

of about 0.1 ohm- cm indicates that Q is larger than conventional 

E 

theory predicts and that is smaller than the base lifetime t_. 

E ^ B 

As discussed in Chapter 5, these new experimental methods will serve 
several purposes in our future research. They help toward determining 
which of the many possible fundamental electronic mechanisms dominate 
in degrading They help enable a better estimate of the maximum 

n achievable from silicon technology. They provide a basis for studies 
aimed to increase and hence ri. 
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The approaches of Chapters 3 and 4 require the validity of 
the shifting approximation of (1) , which implies that a study of 
the dark current-voltage characteristic (V) will suffice in determining the 
origin of the discrepancy. Fundamental ci_-nslU: .tions under- 

lying the validity of (1) are treated in the Appendix. 

Insofar as is possible, each of the Chapters are written as 
independent units to enable the reader with a main interest in a 
later chapter to be able to read that one directly without first 
having to become ■ familiar with the earlier chapters. 
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chapter II. FUNDAMENTAL ELECTRONIC MECHANISMS LIMITING THE 

PERFOSMZVtICE OF SOLAR CELLS 


I. Introduction 

Any attempts to characterize the efficiency of solar cells from a fundamental 
viewpoint must begin with a set of differential equations that describe the 
transport, recombination, and generation of holes and electrons. The efficiency 
obtainable from different material technologies depends upon the fundamental 
limitations imposed on the material parameters appearing in these equations. 

The various basic mechanisms that can determine the properties of these material 
parameters are the main subject of this paper. 

In the sections to follow, two dominant mechanisms are discussed: (a) 

distortion in the energy band, and (b) carrier recombination and generation 
rates or lifetimes. The spatial dependences associated with these mechanisms 
are emphasized — ^both in the direction normal to the surface on which the sun 
shines and that tangential to it. As will be seen, the spatial dependences 
can play a central role in determining the solar cell efficiency. 

II. Basic Mechanisms Affecting Conversion Efficiency 
2.1 The Shockley Equations 

The basic mechanisms that determine the efficiency of a solar cell can 
be analyzed from the basic set of differential equations that govern the transport 
and recombination-generation-trapping of the electrons and holes in a semiconductor. 
These equations were first used by Shockley [1] for thermal recombination and 
generation of electrons and holes. They were generalized by Sah [2] to include 
optical and other generation-recombination processes and further extended to 
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include high recombination-generation or low scattering rates [3] using the 
Boltzmann equation. The complexities introduced by varying material composition 
and heavy impurity doping, resulting in positional dependent electron affinity 
and energy gap, were also discussed and analyzed by Van Overstraeten [4], but. 
•the significance of band and bound states separated by a mobility edge was 
not recognized. The concept of a mobility edge will be further elaborated I 
on in a later section and in a companion paper [5] . For the present purpose 
of delineating the important solar cell parameters, the Shockley equations for 
the specialized condition of d.c. steady-state in compositionally uniform 
material will suffice. These are [2] ; 


J„ = -fqD VN + qy NE (1) 

N n n 

= - qU^NVV^ (la) - 

Jp = - qDpVP + qy^PE (2) 

= - qy PVV (2a) 

P X 

° + ^SS 

° = “ ""^P + ^SS 

V-S = - V^Vp = p/e = (q/e)(P - N + - N^) (5) 


They can be solved for any device geometry and material parameters, as part 
of a boundary-value problem, using the auxiliary conditions 


N = n^exp[q(Vp-Vj^)/kT] , 

(6) 

P = n^exp[q(Vp-Vj)/kT] 

(7) 
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and the boundary conditions to give the total current and applied voltage, 

J - + Jp. \ - Vp(0) - V^(L) (8) 

as a function of such external excitations as be light intensity. In the 
expression for the applied voltage given by (8), .the ohmic contact to the p-type 
is located at x=0 and to the n— type is located at x=L. 

Equation (1) expresses the electron current density, in A/cm^, as the 

sum of a diffusion component, which is proportional to the gradient of the concen- 
tration of the electrons, N, and a drift component, which is proportional to. the 
conductivity, qP^N, and the electric field, E. The coefficient of proportionality 
involves the electron mobility, p^, and diffusivity, D^. Equation (la) gives an 
equivalent way of expressing the electron current density, in terras of the product 
of the electron conductivity, and the gradient of the quasi-Fermi potential 

of electrons, V^. 

Equations (2) and (2a) are the hole current expressions; they are similar 
to those given in Equations (1) and (la) just described for electrons. 

The physical basis for the definition of carrier mobilities and diffusivities 
appearing in Equations (1) to (2a) is a nontrivial one when the band and bound 
states are not separated by a distinct energy band edge. Such a fuzzy band 
edge occurs in heavily doped single crystal and poly crystalline semiconductor 
materials for solar cell applications. A mobility edge is then defined to 
delineate the band states of high carrier mobility from bound states of low or 
negligible carrier mobility. The latter must then be taken into account as 
carrier generation— recombination— trapping centers to be Included in Equations 
(3) and (4). This is further elaborated on in a later section. (See also Sah [5].) 

Equations (3) and (4) are the continuity equations for electron and 
hole currents. They state that the dc steady state is maintained by two components 
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so thac there is no time dependence of electron or hole concentrations. These 
two components are the divergence of the electron or hole current density and 
the net steady-state generation (generation minus recombination) of electrons 
and holes. The latter is written as a volume current source, I„„. It consists 
of all the contributions from the many possible electron— hole generation and 
recombination mechanisms which have been described in detail by Sah [2], For this 
discussion, we shall decompose I into components particularly important for 

M b 

solar cell applications, which can be written as 


"ss ■ 


(9) 


Here, and are the optical generation rate of electrons and holes respectively. 
They are equal, G^ = Gp = G*^, if we consider only the interband optical generation 
of electron-hole pairs from photoexcitation of electrons in the valence band into 
the conduction band and neglect optical excitation of electrons and holes trapped 
at impurity and defect centers. and are the net recombination rates of 
electrons and holes, respectively. 

In solar cells, the term G*^ describes the generation of electron-hole 
pairs .by sunlight. It gives rise to the total photocurrent or short circuit 
current, which is commonly used to characterize the quality of a solar 

Ox.’ 

cell and to guide in the design of efficient solar cells. [6] 

If only the interband optical generation of electron-hole pairs is 
cons.idered as indicated above, then we have the simple steady— state condition 
= Rp = Rgg* The net steady-state recombination rate, Rggj contains 
nonoptical and nontunneling mechanisms discussed by Sah [2]. These mechanisms are 
band and band-bound thermal and Auger-impact transitions which are summarized 
below for the d.c. steady-state condition. 
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( 10 ) 


R, 


ss 


R 


CV 




+ R 


CTV 


The band-to-band (conduction-valence band) net steady-state recombination 
rates are: 


R^^ = q(r^NP - g^) Thermal (11) 

Rpw = (r^ P - g^) + q(r%^N - g%) Auger-Impact (12) 


which can be combined to give 

^CV ~ ^CV ^ ^CV ~ <l(^RR “ s) 

where 


(13) 


■r = r^ + r’^N + r^P, g = g^ + + g^P 


(14) 


Here, r and g are the interband thermal recombination and generation rates 

of electron— hole pairs, r^ is the Auger recombination rate of an electron— hole 

pair in which the energy is carried away by an electron while g" is the impact 

generation rate of an electron-hole pair initiated by an energetic electron. 

r^ and g^ have the same meaning but the third particle is a hole. 

The net' steady-state recombination rate for band-bound or conduction band- 

trap-valence band transitions consists of that for the thermal mechanism, 

T 

Rctv=^SRH considered by Shockley, Read and Hall [ 7 ], and the Auger-Impact 

A 

mechanism, R^.j^* The net d.c. steady-state rate of these two mechanisms for 
single-level generation-recombination-trapping (GRT) center is given by [2] 


R 


CTV 


= R 


CTV 


+ R 


CTV 


qNTT 


c c NP - e e 
n p n p 

c N + e + c_P + e„ 
n n P P 


(15) 
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and the trapped charge concentration, N^, is given by 

c N + e 

T _ n p 

cK+e +cP+e 
TT - n n p p 


The capture and emission coefficients are given by 


( 16 ) 


c 

n 


c*“ + (c”n + c^P) , 
n n 11 


e == + (e^N + e^P) 

n n n n 


(17) 


Cp = Cp + (CpN + c^P), e^ = 6p + (CpN + e^P) 

Here, c^ and are the thermal capture rates and e^ and e^ are the thermal 
n p n p 

emission rates of electrons and holes respectively at a single level GRT center 
whose positional dependent concentration is denoted by N^(x,y,z), c^ is the Auger 
capture rate of an electron by a trap with a second electron carrying away the' 
energy while c^ is the Auger capture rate of an electron by a trap but with a 
hole carrying away the energy. The rates of the corresponding inverse transition 
processes of trapped electrons are denoted by e^ and e^, with the superscripts 
denoting theenergetic electron, n, or hole, p, which cause the impact emission 
transition of a trapped electron. The corresponding rate coefficients for 
Auger capture and impact emission of a hole are denoted by Cp, c^, e^ and e^. 

Equation (5) is the Poisson equation which relates the electric field, E, 
and electrostatic potential, V^, in the solar cell to the macroscopic charge concen- 
trations. The components of the charge concentrations include the concentrations 
of the ionized donors and acceptor impurities, and N^, which may be position- 

dependent, as well as the concentration of the negatively charged GRT centers, 

N^, given by Equation (16). In many of the earlier mathematical analyses of ■ 
semiconductor devices, the trapped charge, -qN^, is neglected. For example, 
Shockley's p-n junction theory [1] did not include the trapped charges in the 


18 



Poisson equation when it was solved to obtain the electrostatic potential in the 
junction, A detailed analysis of the d.c. steady-state trapped charge concentra- 
tion, such as Equation (16), was later given by Shockley and Read in their work 
on recombination kinetics [ 7 ] where detailed analyses of the effect of trapping 
on the steady— state lifetimes were given in the appendix. The importance of 
the trapped charge can be simply illustrated by noting that near the p-n boundary 
of a junction, is nearly zero so that the trapped charge, , can have 

a substantial effect on the electrostatic potential variation in materials with 
non— negligible trap concentration such as solar— grade silicon single crystals 
and polycrystalline films. 

These equations, together with appropriate boundary conditions, 
can determine all aspects of solar cell operation, including solar cell 
energy conversion efficiency. To help achieve better conversion efficiency a 
quantitative understanding is needed of the material parameters appearing in the 
equations, and particularly of the functional dependencies of these parameters. 

In essentially all previous simple analyses of solar cells, as well as in 
detailed numerical computer solutions of these equations for solar cells, 
many idealizing approximations for the parameters have been made, for example, 
the recombination rates R^ and Rp. These approximations were made and then 
propagated to later work for three apparent reasons. First, they simplified 
analysis. Second, the severity of approximation involved has apparently gone 
unrecognized, except in a few instances, and the approximations have been 
treated as if they were-valid. Third, and most important, experimental tools 
have not been developed and applied to supply the needed functional dependencies 
of these parameters in solar cell structures. 
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The functional dependencies of the material parameters therefore govern 
solar cell efficiency. The dependencies have their physical origin in certain 
basic mechanisms as well as in cell inhomogeneities. Among the most important 
of these are: 

(1) Distorted energy bands, 

(2) Interband or band-bound transition or recombination rates, and 

(3) Positional dependences in directions normal and parallel to the 
illuminated surfaces (areal inhomogeneity) . 

We now briefly discuss these factors and their role in determining solar cell 
conversion efficiency for various material systems. 

2.2 Distorted Energy Bands 

At high doping concentrations, and for material with considerable disorder, 
such as polycrystalline films or varying chemical composition, various deviations 
from the idealized band structure can occur, such as: 

(1) The energy band edges are no longer sharp. Electronic energy levels 
can extend beyond the energy positions of band edges of perfect single- 
crystal material, forming thereby so-called band tails [8]; 

(2) The localized impurity states can broaden into an impurity band [9] with 
increasing impurity concentration. At high doping levels,' these 
states can become substantially delocalized; 

(3) A spatial variation of the energy gap can result from macroscopic 
lattice strain Introduced by a high-concentration diffused impurity 
and structural imperfection in the lattice [10] , as well as chemical 
composition variation in graded materials [3,11] and he troj unctions [12]. 
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Each of these deviations tends to produce an effective shrinkage AE^ of 
the energy band gap. Because the intrinsic density n^ depends exponentially 
on the band gap, the gap shrinkage has been imbedded, as a first-order approxi- 
mation, in an effective intrinsic density [13], which then appears directly in 
Eqs. (6) and (7). Note that spatial dependence of the doping concentration 
implies spatial dependence of the effective band gap. 

To examine qualitatively the effect of gap shrinkage on the conduction 
processes described in Eqs. (1) and (2) and on the recombination-generation 
rate in Eqs, (3) and (4), we must distinguish whether the electronic states in 
the presence of high doping are localized or delocalized. Thus, borrowing from 
the concepts used in models describing amorphous semiconductors, we identify 
a critical energy, the Cohen mobility edge [14], which separates localized states 
from delocalized states. The mobility edge plays a dominant role in describing 
charge transport in disordered-aperiodic solids similar to the role the band 
edge plays in periodic solids such as lightly-doped silicon. 

For delocalized states, conduction can occur by drift and diffusion. 

If the variation in doping produces spatial dependence of the band gap and the 
electron affinity, the so-called quasi-electric fields [5,11] produces the drift 
currents of electrons and holes. For localized states, conduction can occur by 
carrier hopping between neighboring states [15] due to thermally activated 
transfer of the trapped charges in the Mott model [16] or by tunneling in the 
Anderson model [17] . 

For ordinary temperatures encountered in solar cell operations, only 
the delocalized or band states will have high carrier mobility and dominate the 
electrical current. In these band states, three effective forces may act to 
produce the flow of electrons and holes. They are: (i) the Coulomb and short- 

range forces, giving rise to the diffusion current due to concentration gradient, 
(ii) the electric force, giving rise to the conventional drift current and (iii) 
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the electric force from the positional dependences of energy gap and electron affin- 
ity, giving rise to additional drift currents. The second and third effective force 
can be combined into an electric field known as the quasi-electric field whose gradient 
gives the drift current components [5] . The total current written in terms of 
the gradient of the quasi-Fermi potentials given by Equations (la) and (2a) 
remain valid as they give zero currents at equilibrium when the quasi— Fermi 
potentials coincide and become spatially constant. 

2.3 Interband and Band-Bound Transitions; Recombination Rates and 

Effective Lifetimes 

The interband and band-bound transition rates of electrons and holes, 
and given in Equation (9), are customarily modeled using oversimplified 
approximations in semiconductor device analysis. For example, they are normally 
written as the excess carrier concentration (lonequilibrium minus equilibrium 
carrier concentration) divided by a constant effective lifetime. This constant 
recombination lifetime model is questionable in the diffused emitter layer of 
single— crystal silicon junction and is likely to be too simplified for polycrys— 
talline solar cell materials in general. For instance, the transition rates or 
the effective lifetimes of electrons and holes depend on the several energy- 
momentum conservation mechanisms such as the phonon-assisted or thermal, 
the photon— assisted or optical and the three— body or Auger— Impact processes, [2] 

The recombination rates of such transitions depend on the density, the 
energy levels and the electron and hole capture and emission rates of the 
localized centers , and on the positional dependences of the energy 
band gap. These rates have already been summarized in Equations (10) to (17), 

Only under rather special circumstances, namely low illumination and extremely 
high illumination intensities, can these rates be characterized by constant 
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lifetime provided the concentrations ot ,the recombination centers' are 
spatially constant [6]. The lack of spatial variations is most unlikely in 
polycrystailine solar cells but also in single crystal cells in view of its 
large geometry as well as the heavily disordered diffused layer. A major 
reason for the use of the constant lifetime approximation in the past was the 
unavailability of accurate experimental data on the recombination parameters 
and the lack of sensitive experimental techniques which can determine these 
parameters in solar cells instead of the uniform semiconductor bulk. New 
methods have now been developed so that one can measure these parameters in 
solar cell structures [6, 18, IS]. 

Aside from the mechanisms mentioned, another type exists that can influence 
the effective interband transition rates. This involves the tunneling transition 
via deep'-level defect centers [2], Among the many transition processes of this 
general type [20] , we shall describe here a particular two-step process. This 
process assumes that an electron in the conduction band is first captured by a 
localized deep level and subsequently tunnels to the valence band or to the 
valence-band tail states. It has been demonstrated experimentally that this is 
the excess-current mechanism in gold-doped silicon tunnel diodes [21 ] . 

In polycrystailine films, recombination processes occur at grain boundaries 
by any of the energy /momentum conservation mechanisms mentioned previously. 

As a phenomenological first-order approximation, one can characterize this 
recombination in terms of an effective surface recombination velocity S acting 
over the surface area of the grain boundaries. Choice of various assumed 
structures of the grains can then enable theoretical study, for example, of the 
relation between the size of the grain and the effective lifetime of the 

t 

polycrystailine material. [22] 
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opaLXcix uepenoences 


The effect of the mechanisms just discussed on the solar-cell efficiency 
depends critically on the spatial dependences of these mechanisms in the direction 
normal to the surface of the cell, as well as in the transverse direction. These 
are discussed in the following two sections. 

Variations Normal to the Cell Surface 

The shrinkage of the energy gap may become large in regions that are 
heavily doped, for example, in the diffused emitter region of a pn junction 
solar cell. There it acts to increase the minority charge stored. The 
spatially variation impurity concentration in the diffused emitter will 
also produce a built-in electric field which can reduce the transit time 
required for the minority carriers to cross the emitter region. Both of 
these effects tend to increase the dark saturation current, thus lowering 
the open-circuit voltage and the conversion efficiency. 

In the heavily-doped emitter of a p-n junction solar cell, one expects the 
defect density to rise sharply near the degenerately-doped surface. This 
expectation comes from the experimental evidence [23] on single-crystal silicon 
which showed that during silicon' device processing, recombination centers are 
generated at the strained and disordered surface layer and diffuse rapidly into 
the bulk. 

These experiments further suggest that the centers are silicon vacancies 
or vacancy complexes. Thus, using the analogy to chemical reactions reviewed 
by Fuller [24], one finds the following first-order model for the dependence of 
the recombination center density on the concentration .of the impurity 
in the diffused emitter and space-charge regions: 

Ntt(x) = K[Nj^jj(x) + N^] . 

In developing fhis model, one assumes approximate thermal equilibrium during 

the diffusion of impurities. This justifies the law of mass action, which one 

combines with the electrical neutrality condition. The different values of r are 
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derived from different hypotheses about the dominant chemical reaction. For 
reactions corresponding to the formation of divacancies, the choice r = 2 
provides an approximate fit. Experimental techniques exist [6] by which one 
can study the validity of the theoretical anticipation given in Equation (18), 
and determine the value of the parameter r if the theoretical anticipation is 
valid for the cell under study. 

The occurences of spatial dependences and high concentrations of 
recombination centers would degrade solar cell performance in various ways. 

High recombination rate in the junction space-charge region and in quasi- 
neutral emitter layer of p-n junction solar cells increase the dark saturation 
current and thereby lower the open-circuit voltage. Moreover, recombina- 
tion in the space-charge region, where electron and hole concentrations are 
large, tends to decrease the fill factor as well as the short circuit current 
[6]. Thus, high recombination rates in the junction space-charge region and 
in the quasi-neutral emitter act in three ways to decrease the conversion efficiency. 
In addition, if the lifetime in the quasi-neutral base were to decrease with 
position in the direction toward the surface, then the dark saturation current 
would further increase. Hence, this too can contribute to lowering the 
conversion efficiency. 

Variation Across Surface (Areal Inhomogeneity) 

We have thus far considered a one-dimensional model of the cell in which 
the only coordinate of interest has been that measuring the distance from the 
surface. But the solar cell is a large-area device, and inhoraogeneities across 
its area could play a significant role in governing the performance. In 
particular, we note the existence of a spatial distribution of impurity clusters, 
thermodynamically stable, occurring in heavily-doped emitter regions [25] • Thus, 
the inhomogeneities may result, in part, from inescapable statistical 
fluctuations and, in part, from man-made contributions that could be minimized. 
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One can systematically view the influence of the one-dimensional mechanisms 
discussed previously in the context of inhomogeneity across the area of the cell. 
We stress that the overall effect on solar cell characteristics is not -a simple 
average over the area of the cell. Indeed, a small defective area could dominate 
in determining the conversion efficiency of a large area solar cell. For example, 
certain of the mechanisms are much emphasized in a region of large-impurity 
concentration due to clustering. An example is energy band gap shrinkage, 
which enters as an exponent in determining cell properties in a localized region; 
another example is the associated increased defect density and the increased 
recombination rates corresponding to it. The model of a solar cell as a 
collection of mini-cells connected in parallel with one another, although 
much simplified, qualitatively suggests the trends in behavior introduced, by 
areal inhomogeneity-; that those mini-cells with relatively high doping or 
defect density (or with small grain size in a polycrystalline film) can severely 
degrade the performance of the device. Hence, areal inhomogeneity could play 
a dominant role and establish a technological limitation on the solar conversion 
efficiency obtainable. 

Ill . SUMMARY 

Certain basic mechanisms in semiconductor material and their dependence 
on position were the subjects of this paper. In outline form, the subjects 
treated were: 

(a) Distorted and spatially varying energy band in highly-doped and variable 
chemical composition material 

(b) Bulk, surface and grain boundary recombination 

(c) Spatial dependencies normal to the surface and areal inhomogeneities 
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Each of the above can strongly influence the performance of solar cells. 

For example, energy band distortion present’ in a highly-doped emitter can limit 
the open-circuit voltage obtainable in p-n junction solar cells. High recombina- 
tion -rates can affect the short-circuit current, the fill factor, and the open- 
circuit voltage. The spatial dependence of the energy band distortion and the 

I 

recombination rates, both across the area of the cell and in the direction normal 
to the- surface, can determine the degree to which the fill factor, the short- 
circuit current, and the open-circuit voltage are degraded. Recently, experimental 
methods have been developed and applied to silicon solar cells to determine the 
species of recombination centers [6] and the energy gap shrinkage and the 
effective lifetime in the heavily-doped emitter [19]. 

All of the mechanisms above apply both to single-crystal as well as to 
polycrystalline thin-film materials. For polycrystalline cells, processes at 
grain boundaries can play a central role in determining the solar conversion 
efficiency [22]. 

Research studies are needed to improve our understanding of the basic 
mechanisms treated here. Detailed quantitative studies are required to assess 
the ultimate capabilities of any particular material technology. Thus such 
research can potentially yield a basis for choosing from among alternative 
technologies. Moreover, for solar cells made from any given technology, such 
as the technology of single-crystal silicon, studies of the basic mechanisms 
can help guide the fabrication processing toward maximizing the conversion 
efficiency. 

In chapters 3 and 4, we describe studies of silicon devices -that help 
toward achieving all of these objectives. 


27 



REFERENCE FOR CHAPTER 2 


[1] W. Shockley, "The theory of p-n junctions in semiconductors and p-n 
junction transistors," Bell Syst» Tech. J . , vol. 28, pp. 435-489, 1949. 

[2] C. T. Sah, "Equivalent circuit models in semiconductor transport for 
thermal, optical. Auger-impact, and tunneling recombination-generation- 
crapping processes," physica status solidi a, vol. 7, pp. 541-559, 1971. 

[3] C. T. Sah and F. A. Llndholm, "Transport in semiconductors with low 
scattering rate and at high frequencies," Solid-State Electronics , 
vol. 16, pp. 1447-1449, 1973. 

[4] R. J. Overstraeden, H. J. DeMan and R. F. Mertens, "Transport equations 
in heavy doped silicon," IEEE Trans. Electron Devices , vol. ED-20, 

pp. 290-298, 1973. 

[5} C. T. Sah, "Carrier generation, recombination, trapping and transport in 
semiconductors with position dependent composition," IEEE Trans. Electron 
Devices , this issue. 

[6] C. T. Sah, "Detection of recombination centers in solar cells from 
• capacitance transients," IEEE Trans. Electron Devices, this issue. 

[7] W. Shockley and W. T. Read, Jr., "Statistics of the recombinations of 
■ holes and electrons," Phys. Rev ., vol. .87, pp. 835-842, 1952. 

[8] I. M. Lifshitz, J. Exp. Theor. Phys ., vol. 2, pp. 117, 137, 156, 1942. 

[9] _ G. W. Castellan and F. Seitz, "The transition from insulating to metallic 

behavior in semiconducting silicon," Phys , Rev . , vol. 79, pp. 216-222, 
1950; C. S. Hung, "Theory of resistivity and Hall effect at very low 
temperature," Phys. Rev ., vol. 79, pp. 727-728, 1950; C. Erglnsoy, "On 
the mechanism of impurity band conduction in semiconductors," Phys. Rev. , 
vol. 80, pp. 1104-1108, 1950; and T. N. Morgan, "Broadening of Impurity 


28 



bands in heavily doped semiconductors," Phys. Rev ., vol. 139A., pp. A343- 
A348, 1965. 

[10] V. Fistul, Heavily Doped Semiconductors , Plenum Press, N.Y., 1969. 

[11] H. Kroemer, "Band structure of semiconductor alloys with locally varying 
composition," Bull. Amer. Phys. Society , vol. 1, p. 143, Paper V9, 1956; 
"Quasi-electric and quasi-magnetic fields in nonuniform semiconductors," 

RCA Review , vol. 28, pp. 332-342, 1957. 

[12] A. G. Milnes and D. L. Feucht, Heteroj unctions and Metal-Semiconductor 
Junctions , Academic Press, 1972. 

[13] D. D. Kleppinger and F. A. Lindholm, "Impurity concentration dependence 
of the density of states in semiconductors,” Solid-State Electronics, 

vol. 14, pp. 199-206, 1971; D, D. Kleppinger and F. A. Lindholm, "Impurity • 
concentration dependent density of states and resulting Fermi level for 
silicon," Solid-State Electronics , vol. 14, pp. 407-416, 1971; H. J. DeMan 
"The influence of heavy doping on the emitter efficiency of a bipolar 
transistor," IEEE , vol. ED-13, pp. 833-834, 1971; R. J. Van Overstraeten, 

H. J. DeMan, and R. P. Mertens, "Transport equations in heavy doped silicon," 
IEEE , vol. ED-20, p..-290, 1973; H. J. DeMan, R. P. Mertens and R. J. 

Van Overstraeten, "Influence of heavy doping effects on the f;p prediction 
of transistors," Electronic Letters , vol. 9, pp. 248-249, 1973; R. P. Mertens 
H. J. DeMan and R. J. Van Overstraeten, "Transport equations in heavy doped 
silicon," IEEE Trans. Electron Devices , vol. ED-20, pp. 772-778, 1973; 

M. S. Mock, "Transport equations in heavily doped silicon, and the current 
gain of a bipolar transistor," Solid-State Electronics , vol. 16, pp. 1251- 
1259, 1973; M. S. Mock, "On heavy doping effects and the injection efficiency 
of silicon transistors," Solid-State Electronics , vol. 17, pp. 819-824, 1974; 
and M. P. Godlewski, H. W. Brandhorst, Jr., and C. R. Baraona, "Effects 
of high doping levels on silicon solar cell performance," High Efficiency 


29 



silicon Solar Cell Meeting , NASA-Lewis Research Center, 1974; Record of 11th 
Photovoltaic Specialists Conf ., pp. 32-36, 1975. 

[14] M. H. Cohen, H. Fritzsche, and S. R, Ovshinsky, "Simple band model for 
amorphous semiconducting alloys," Phys. Rev. Lett ., vol. 22, pp. 1065- 
1068, 1969. 

[15] E. M. Conwell, "Impurity band conduction in germanium and silicon," 

Proc. IRE , vol. 103, pp. 51-61, 1956. 

[16] N. F. Mott, "On the transition to metallic conduction in semiconductors," 
Can. J. of Physics , vol. 34, pp. 1356-1368, 1956. 

[17] P. W. Anderson, "Absence of diffusion in certain random lattices," 

Phys. Rev , vol. 109, pp. 1592-1605, 1958. 

[18] ■ C. 1. Sah, L. Forbes, L. L. Rosier, and A, F, Tasch, "Thermal optical 

emission and capture rates and cross sections of electrons and holes at 
imperfection centers in semiconductors from photo and dark junction 
current and .capacitance experiments," Solid-State Electronics , vol. 13, 
pp. 759-788, 1970; For a review of the applications of these techniques to 
recombination centers in silicon, see C. T. Sah, "Bulk and interface 
imperfections in semiconductors," Solid-State Electronics, vol. 20, pp. 1-20. 
1976. See also reference 5. 

[19] (a) F. A. Lindholm, A. Neugroschel, C. T. Sah, M. P. Godlewski, and 

H. W. Brandhorst, Jr., "Methodology for the experimental determination 
of gap shrinkage and lifetimes in the emitter and base of pn-junction 
solar cells," IEEE Trans. Electron Devices , this issue; also Record of 
12th Photovoltaic Specialists Conf ., Nov. 1976, 

(b) F. A. Lindholm, A. Neugroschel, and C. T. Sah, "Junction modeling 

for solar cells — theory and experiment," Digest of 1976 Int~. Electron 
Devices Meeting , Dec. 1976. 


30 



(c) A. Neugroschel, F. A. Lindholm, and C. T. Sah, "A method for 

determining the emitter and base lifetimes in pn-j unction diodes," 

IEEE Trans. Electron Devices, submitted for publication. 

[20] C. T. Sah, "Excess current in semiconductor tunneling," in Tunneling Phenom- 
ena in Solids, Ch, 14, (Editors: E. Burstein and S. Lindquist) Plenum Press, 
N.Y., 1969. 

[21] C. T. Sah, "Electronic processes and excess currents in gold— doped narrow 
silicon junctions," Phys. Rev, pp. 1594-1612, 1961. 

[22] C. T. Sah and F. A. Lindholm, "Characteristics of solar cells on granular 
semiconductors," Proceedings of the 12th IEEE Photovoltaic Specialists 
Conference , Nov. 15, 1976. 

[23] C. T. Sah and C. T. Wang, "Experiments on the origin of process induced 
recombination centers in silicon," J. Appl. Phys ., vol. 46, pp. 1767-1776, 

’ 1975. 

[24] C.'S. Fuller in Semiconduc tors , Ch. 5, (Editor: N. B. Hannay) Reinhold Corp., 
N.Y., 1959. 

[25] W. Shockley, "Problems related to p-n junctions in silicon," Solid-State 
Electronics, vol. 2, pp. 35-67, 1961. 


31 



CHAPTER III. A. METHOD FOR DETEEMINING THE EMITTER AND BASE' LIFETIMES 


I . INTRODUCTION 

Past work on p-n- junction diodes has not provided a means -for 
the experimental determination of the emitter lifetime. This has 
prevented a quantitative assessment, based on experiment, of the 
importance of carrier recombination in the emitter to the electrical 
characteristics of p-n diodes, solar cells, and bipolar transistors. 

A main purpose of this paper is to emphasize and to demonstrate, 
by Q.xpejLOnzyvt , the significance of the heavily-doped emitter region of a diode 
to its forward current- voltage characteristics. The experimental evidence 

to be presented shows that the emitter can play a significant role in 
the response to both static and transient excitations. In -the course 
of the study, the minority- carrier lifetimes are found experimentally 
for the emitter and the base regions. As will be seen, the determina- 
tion of these two lifetimes demands that more care be given to 
constructing experiments and interpreting experimental evidence than 
has been given previously. A second main purpose of this paper is to 
describe the techniques by which we distinguish the emitter lifetime 
from the base lifetime. 

The reference to emitter and base used above for diodes borrows 
the terminology used for transistors. In the paper, this terminology 
will refer to a p-n diode consisting of a heavily-doped p-type emitter 
and an n-type base that is more lightly doped. The emitter and base 
are quasi-neutral regions separated from each other by the space-charge 
region of the jtmction'. For concreteness the treatment throughout this 
paper will refer to a p^-n structure; an analogous treatment applies 
for n^-p structures. 
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The significance of' the emitter applies not only to diode theory but 
also to the theory of operation of other devices, such as the p-n junction 
solar cell and the bipolar transistor, whose structure involves one or more 
p-n junctions. For example, the physical processes occurring in the emitter 
of a bipolar transistor c^ limit the achievable maximum common-emitter 
current gain. For the p-n junction solar cell, as Chapters 1 and 2 have 
discussed, these processes can limit the achievable maximum solar energy 
conversion efficiency. The experiments and particularly the methods for 
interpreting experimental data treated in this paper for diodes, therefore, 
apply also to junction transistors and solar cells. A key connection between 
p-n-junction diodes and p-n-junction solar cells is the shifting approxima- 
tion (Eq. 1, Ch. 1), the validity of which forms the subject of the Appendix. 

The emitter will be significant to the behavior of a device whenever 
the current that supports the recombination in the emitter constitutes an 
appreciable fraction of the currents flowing at the device terminals. A 
large emitter recombination current arises from high recombination rates in 
the emitter. As has been discussed in Chapter 2, such rates can be traced 
to various physical origins, which fall into two broad categories [11: 

(a) ban.d~zd.gQ. dlMtofvtion that may be present in heavily-doped 
material can affect recombination, diffusion and drift; and 

(b) crystalline imperfections, formed near the emitter surface 

and diffused into the emitter during high- temperature 
fabrication, can act as H.Z(lombZnatAMn producing 

a lower lifetime in the emitter than in the base. 

The effect of recombination in the emitter, category (b) above, receives 
the emphasis in this paper. The design of the diode structures used here 
purposefully makes negligible the effects of band-edge distortion. Chapter 
4 will treat the effects of band-edge distortion, category (a) above. 
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Section II gives details about the design and fabrication of the diodes. 
Section III gives the values of the recombination-generation parameters of 
the deep energy levels present in the base, as revealed by transient- 
capacitance experiments. The energy levels studied in detail fall in the 
upper half of the bandgap. These energy levels control the lifetime of 
minority carriers in p-type material. Thus studies of the type reported 
here, made on p-on-n diodes, can potentially yield the parameters controlling 
and limiting the base lifetime in n-on-p solar cells . This relationship is 
noted because n-on-p solar cells show advantages for space applications, 
and are the type of cell of main interest in this NASA grant research. 

Section IV characterizes the static current-voltage dependence in a way 
that is consistent with the underlying physics. Such a characterization is 
essential to 'the proper calculation of values of lifetimes and of other 
material properties. Section V demonstrates and discusses erroneous con- 
clusions about the base lifetime calculated from improper current-voltage 
characterizations used in the past. In Section VI, methods for interpreting 
data are described that distinguish the material properties of the emitter 
from those of the base. Applying these methods then yields the values for 
the lifetimes of the emitter and the base, which are tabulated and discussed 
in Sections VI and VII. 

II. DIODE FABRICATION AND PHYSICAL MAKE-UP 
Boron was diffused into n-type silicon substrates of (111) orientation 
to form p -n junctions. Four sets of diodes were made, each with a different 
phosphorous concentration N^^^^ in the substrate, ranging from 10^^ to 10^^ cm 
Various substrate doping concentrations were used to permit the testing 
of qualitative expectations about the significance of the emitter. On quali- 
tative grounds, one notes, for example, that the excess minority charge 
stored' in the base is inversely proportional to -the doping concentration of 
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the base (substrate) , as a first approximation. Hence the higher the 
substrate doping the more one might expect the emitter to dominate the 
current- voltage characteristics of the diode. This expectation is borne 
out by the experimental results described in Sections V and "VI. 

In the fabrication of the diodes, after standard surface cleaning 

[2], a 0.5 ym thick oxide layer was grown in a wet oxygen ambient. Holes 

2 

of rectangular area 25x42 mil were opened in the oxide on the top sur- 
face using standard photolitographic techniques. Boron was then pre- 
deposited into these oxide holes from a solid planar Boron-plus source 
■at 1100®C for 30 minutes in mixed + 0 ^ carrier gas. The boron glass 
was then etched away and a boron drive-in diffusion was done at 1200°C 
for three hours in dry oxygen. 

Slow pulls lasting about 5 min. were used during the predeposition 
step, mainly to avoid thermal warpage of the diffusion sources. The pull 
rate was high during the drive-in step. For the studies of interest here, 
pull rates are not very important since only extremely slow cooling could 
reduce the recombination-center densities in the bulk [3] . 

Contact holes to the p"^ region were then opened in the 0.35 pm oxide 
layer grown during the drive-in step. Aluminvim was evaporated on the 
top side and was selectively removed using photolitrography to define 
contacts to the p region and to create an MOS guard-ring gate overlapping 
the edge of the diffused layer. Aluminum contacts were sintered at 450° 
for 30 minutes in dry N^. The back side was lapped to assure a good ohmic 
contact to the substrate. The finished wafers were then- scribed into 
125 mil- square chips and bonded to TO-5 headers at 400°C. Gold wires 
were, bonded to the aluminimi contact pads and- to the MOS guard-ring. 

The phosphorus concentration in the substrate was obtained with good 

agreement from two standard methods: from the dependence of capacitance 

2 

(I'/C ) on reverse voltage and from four-point probe measurements of the 
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starting ras\ sLi vit^ . VaUu's obta.iiu'vl liy Uu' ca|viei t .iiuu' inrl lu'vl u'ri'p 
then used because they are more accurate. I3y measurement oi' sheet resist- 
ance using a four-point probe, the boron surface concentration was found 
19 -3 

to be about 4x10 cm . The j\mction depth was determined by angle lapping 

and staining. The junctions are deep; for each base doping concentration, the 

emitter thickness exceeds seven microns. This part of the design 

is meant to reduce the effects of band-edge distortion to secondary 

importance, as will be explained in Section 6.1. 

Table I summarizes the physical make-up of the devices. The values 

of the Fermi level E_ given there are calculated from 

1 ? 

Ep - Ej = kT/q ln(Nj^j^/n^) 

10 -3 

using n^ = 1.08x10 cm at 297. 5°K. 

III. measurement of deep levels in the base 

To determine the extent to which the emitter influences the per- 
formance of the diodes just described, the contribution from the base 
must be understood and determined. The component of current coming 
from the base arises from the recombination of holes with electrons 
through deep- level recombination centers in the silicon energy gap. 

The properties of these centers are described by the following para- 
meters: N , the concentration of the centers, e and e , the elec- 

iJ- n p 

tron and hole emission rates, c and c , the electron and hole cap- 

n p ^ 

ture rates. For low-level injection, recombination is characterized 

by the minority-carrier lifetimes, t = 1/c N„^ for p-type material 

no n TT c- j jr 

^po ~ ^^'^p^TT n-type material. Various methods [4,5] exist 

for measuring the parameters, e , e , c , and c , including their 

TT n p n p 

positional dependence. 
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TABLE I 


Summary of Devices Studied 


Device No. 

Base doping 

NDD(cm-3) 

Junction depth 
Xj (ym) 

E -E (297. 5»K) 
(eV) 

2-1-lpl 

1.5 xlO^^ 

10.8 

0.245 

2-2-lpl 

1.34x10^^ 

9.6 

0.300 

2-4-lpl 

1.12x10^^ 

8.5 

0.355 

2-5-lp5 

6.7 xlO^^ 

7.5 

0.401 
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•For our diodes, v;e studied the energy levels in the band gap by the 
capacitance- transient method [4,5]. Use of this method indicated that 
there are no levels in the lower half of the gap. In the upper half, 
the same four levels were observed, independent of the substrate phosphorous 
concentration. The thermal activation energies, densities 

N^T were determined to be as follows: 


93 

meV 

= 3.3x10^^ 

-3 

CII} 

270 

meV 

12 

4.3x10 

-3 

cm 

285 

meV 

2.0x10^^ 

-3 

cm 

542 

meV 

3.5x10^^ 

-3 

cm 


The thermal activation energies were obtained by fitting the experimental 
temperature dependencies of the thermal emission- rates to the Arrhenius 
equation, as in Appendix A. 

The electron thermal- emission rates for all four levels, determined 
by the transient capacitance method, are listed in Appendix A toge^er 
with the hole thermal-emission rates for the midgap level, as calculated 
from the temperature dependence of the reverse-biased current [4] . 

The energy levels at 93 and 285 meV are probably due to the inter- 
sticial Mg double donor center [6] , since the boron diffusion source 
contains MgO as one of its main constituents. The impurities responsible 
for the other two centers were not identified, but they also probably 
come from the binder impurities in the boron diffusion source. The mid- 
gap level at 542 meV has the largest density and will thus tend to 
dominate in determining the base recombination current. In all four 
devices, the density of each level listed above was approximately 
the same, within a factor of two. Measurements also indicated that 
was independent of position in the base. 
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Although it is a well-3cnown procedure [4] , we note again that the 

MOS guard-ring has to be a.c. grounded during capacitance transient 

* 1 * 

measurements to avoid capacitance coupling between the p and n regions 

through a floating MOS aluminum gate on the surface. Since- the 

capacitance measurement is a three-terminal measurement, the MOS 

capacitance to ground is not recorded. Note also that levels in the 

lower half of the band gap will give a decreasing capacitance transient 
t t ' 

because e > e , which contrasts with the increasing capacitance 
p n 

transient due to the levels in the upper half of the band gap, which 

t t 
have e > e . 
n p 

During measurements each device was mounted on a cold finger 
inserted in a liquid nitrogen dewar. Constant temperature was maintained 
by a temperature controller of a type designed by Sah and co-workers [7] . 
The temperature was calibrated to an accuracy of ±0.2°C, and the tem- 
perature change during the measurement was less than ±0.1°C. 

IV. CURRENT- VOLTAGE CHARACTERIZATION 

The determination of material parameters from measured data requires 
that the data be treated and the behavior of the device be described in 
ways that are consistent' with the physics underlying device operation. 

The physical mechanisms governing the operation of a forward-biased 
junction diode produce a current consisting of three components: 

(a) A diffusion component due to the diffusion and recombination 
in the two quasi-neutral regions of the junction, as obtained 
. by Shockley in the ideal diode theory [8] , and extended by 
Moll and Ross [9] to include the effects of drift due to a 
built-in electric field. 
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(b) A recombination component: due to carrier recombination in 
the bulk space-charge layer of the junction [10] and 

(c) Surface components , including the surface channel current, 
the recombination current through bulk centers in the sur- 
face space-charge layer as well as through surface and*, 
interface states [11] . 

In silicon diodes, the recombination current in the bulk space- 
charge layer (component b above) dominates at small forward bias for 
room temperature or below while the diffusion current (component a 
above) dominates at medium bias levels [10] . Surface currents (components 
c above) obscure the properties of interest to this study, and the MOS 
guard-ring gate of the diode structure described in Section II provides 
a .means to diminish and nearly eliminate their contributions. The mag- 
nitudes of the surface components depend on the surface potential, which 
can be controlled by applying a voltage to the MOS guard-ring gate elec- 
trode [11] . By adjusting this voltage to minimize the diode current, 
the influence of the surface components can be suppressed and in some 
devices made negligible compared with the bulk components (a) and (b) above. 

The total diode current is the sum of all three components. The 
recombination and surface components (components b and c) can be com- 
bined into one term, l^ = [exp(qV/m^T) -1], which is an empirical 
relationship that applies to the devices studied here. The total current 
is then the sum of two terms, given by 

I = I^^[exp(qV/m^kT) - 1] + I^^^[exp(qV/kT) - 1] (1) 

where and m^ are the coefficient and the reciprocal slope of the 
surface and bulk space- charge- layer components (components b and c) . 
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In Eq. (!) , the second term is the diffusion current (component a) for low-level' 


injection conditions. The coefficient which is customarily 

called the saturation current, consists of two parts, one from carrier 
recombination, diffusion and drift occurring in the quasi-neutral base 
and the other from these same processes occurring in the quasi— neutral 
emitter . 

The expression for the diode current given in Eq. (1) correctly describes’ 
the physics underlying device behavior for low- level- injection conditions. 
Thus, as we shall see, experimental study of the saturation current 
can help disclose the values of the material parameters of the quasi- 
neutral regions, including the values of the emitter and base lifetimes. 


V. INCORKECT PROCEDURES FOR DETERMINING LIFETIME 

Any characterization using an expression different from Eq. (1) is 
incorrectly linked to the underlying physics for low-level conditions and 
can lead to erroneous conclusions about the material parameters of the 
diode. As an example of such an erroneous link to the physics, consider 
the characterization [12,13]: 

qV/m kT qV/m kT 

I = ^01^® - 1) + - 1) f 1 < m^ < m^ (2) 

The first teinn in this characterization is an empirical fit which has 
been regarded as related to the diffusion current (component a) ; in this 
empirical fit frequently is greater than 1. The second terra is an 
empirical fit regarded as related to the bulk and surface space-charge- 
layer currents (components b and c) . In some past work [14] , the 
reciprocal slope m^ has been set equal to 2, which corresponds to a special 
case of assumed properties of the recombination centers. 
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Equation (2) can describe the terminal characteristics of the single 
diode for which the empirical fit is made. But if this expression is 
used to calculate such material parameters as the base lifetime, it can 
yield values that differ by orders of magnitude from the actual values. 

To illustrate this point, consider the forward-current characteristics 
shown in Fig. 1 for the diodes described in Section II. The characteristics 
were measured at a temperatxire of 297. 5°K, and the surface component of 
the current was minimized by adjusting the voltage on the MOS guard-ring 
gate. The logarithmic dependence of current on voltage shows non-unity 
slope for all four types of diodes shown in Fig. 1, Most silicon diodes 
show a non-unity slope. 

In an incorrect procedure of evaluation used by some authors [12,13], 
the diffusion current is approximated by the first term in Eq. (2) . Then 
the saturation current I^^ is the intercept on the Jlnl axis obtained by 
extrapolating the line of reciprocal slope m^. The values of 1^^ resulting 
from this procedure, together with the reciprocal slope m^, are listed in 
Table II. In the second step in the incorrect procedure, the current 1^^ 
is related to the material parameters by the Shockley ideal formula [8] , 


Iqi = Aqn. 


/S 




+ 




BOl 


N ' T 
AA EOl-J 


(3) 


where A is the device area, n^ is the intrinsic carrier concentration in 
silicon, Dg, are doping density, diffusion coefficient and 

minority carrier lifetime, respectively, in the base, and N^, D^, 
are the same parameters for the emitter region. 

The incorrect procedure continues by assuming that only the base 
is important to the observed behavior, which gives 
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TABLE II 




Summary of Material 

Parameters 

at 295. 5“K 




Device No. 

A 

^01 

A) 

B 

”l 

C 

(nsec) 

D 

T (JCR) 
m 

(nsec) 

E 

T (OCVD) 
m 

(nsec) 

F 

"^B 

(nsec) 

G 

"^B 

T (OCVD) 
m 

2-1-lpl 

500 

1.11 

4.3 

255 

325 

330 

“1.0 

2-2-lpl 

52 

1.09 

3.8 

150 

115 

120 

“1.0 

2-4- Ipl 

11 

1.14 

0.9 

90 

115 

120 

“1.0 

2-5- lp5 

4.6 

1.15 

0.08 

40 

105 

i45 

1.4 






( 4 ) 



The incorrect base lifetime is then solved from Eqn. (4) for each- 

yielding the values in Table II listed in Column C as Note 

that the values deduced by this incorrect procedure are unreasonably 
small, especially for the more heavily doped devices. As will be seen 
shortly, the key to the correct procedure and the correct link to the 
physics lies in the realization that, in accord with Eq. (1) , the line 
used to determine the saturation current must have unity slope. Extra- 
polation of this line will then give much larger values for the base 
lifetime t than -those listed in Col\amn C of Table II. The correct 
procedure of requiring unity slope for the ideal Shockley component of 
current has been known for nearly twenty years. 

A further discrepancy resulting from -this incorrect procedure is 
disclosed by comparing the values of in Table II, found from the 

I-V curve, with values of determined by more direct measurements of 
the lifetime. There are several standard methods of direct measurement 


of the base lifetime in p-n- junction diodes. The methods used most 
often are the junction- current- recovery (JCR) [15] and open-circuit- 
voltage-decay (OCVD) experiments [16] . For both methods, the conventional 
assumption is made -that most of the excess charge resides in the base 
region, and, hence, the result of these measurements is generally 
assumed to give the effective base lifetime. In Table II, we compare 


with T (JCR) and x (OCVD). Substantial disagreement exists. The 
BOl m m 

smallest value of the ratio x /x„_, is about 30. 

m BOl 


The incorrect procedure [12,13] just described is based on the 


following two assumptions; 
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(a) The saturation current can be determined by extra- 

polating the measured dependence of Jlnl versus qV/kT using 
the empirical fit of Eq. (2) , with slope > 1. 

(b) The contribution to diode current from carrier recombina- 
tion, drift and diffusion in the emitter is negligible 
compared with the contribution of these processes in the 
base. 

Assumption (a) is incorrect for all p-n- junction diodes. It is not 
correctly based on the device physics. Assumption (b) is incorrect 
For diodes with heavily-doped bases. Its use can introduce sizable 
error in calculating the base lifetime and other device material 
parameters, as we shall see. 

VI. CONTRIBUTION OF THE EMITTER REGION AND DETERMINATION OF BASE AND 

EMITTER LIFETIME 

To determine the contribution of the emitter requires a model that 
includes the recombination, diffusion and drift processes which occur 
in the quasi-neutral emitter region. As Eq. (3) indicates, the original 
diode theory of Shockley [8] accounted for the recombination and dif- 
fusion of carriers in the emitter. But the Shockley model did not ex- 
plici tly account for several mechanisms that may be present in modern 
silicon diodes. These mechanisms arise from various sources. For 
example [1] : 

(a) the large electric field due to the gradient of the impurity 
concentrations produces minority-carrier drift as well as 
diffusion, which may be thought of as enhanced diffusion; 
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(b) the minority- carrier lifetime may depend on position, 
dccrpas.i ng sharj>ly in the direction toward the 
<iocjunoratcly-dof)od om.ittor surFaco duo to an increasing 
density of recombination centers toward the surface; 

(c) the distortion of tlie energy band edges in the highly- 
doped material' may also influence the spatial dependence 

i 

of the carrier recombination rates in the emitter. 

We now outline a simple lumped model [17] that is general enough to 
include all of these mechanisms, and discuss the use of this model 
in determining device material parameters, including the emitter and 
base lifetimes. 


6 . 1 Model Including the Effect of the Emitter 

As in Eq. (1) , the current in a diode consists of components from 
the qua si -neutral emitter and base regions together with components 
coming from the bulk and surface space-charge regions. Using the 
principles of charge control [18] , we rewrite Eq. (1) as 

■ ^E'^'^E ^ ^^sc'^sc^ . 

in which and are the excess- minority-carrier- charges in the quasi-neutra 
emitter and base, Q represents the effective excess mobile 
charges of the bulk space-charge region and the surface region, and 
T„, To' Trv- respective charge-control time constants or 

relaxation times. For time-varying conditions, the quasi-static ap- 
proximation customarily used in charge-control analysis gives 


% 




i (t), = — + — + 

dt 

E B 




dt 




^SC 

,T 

^ SC 




sc 


dt 


( 6 ) 
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for the diode current- ■ To simplify the. model, we neglect the last 
two components, which are bracketed in. Eq, (6) . As has been discussed, 
the surface currents can be suppressed by applying a voltage to an MOS 
guard-ring. The current from the bulk space-charge region and the rem- 
nants of the surface current not disposed of by the guard-ring can be 
removed from the experimental data by methods to be described. 

If the diode current arises only from proce.sses occurring in the 
quasi-neutral regions of the base and the emitter, then the .expi^ession 
for the static current given in Eq.' (1) reduces to 

I = I^^Q[exp(qV/kT) - ’I] , (-7) 

and the excess minority carrier charges, Q and Q , in the quasi-neutral 

B 

emitter and base depend on the diode voltage V according to' 

^E " - 1)]. Qg = QgQ[exp(qV/kT) - 1)] (8) 

The validity of Eqs. (7) and (S) require the additional assumptions that 
no carrier high- injection levels occur in the quasi-neutral regions and 
that the effect of series resistance is negligible. These as.sximptions 
also necessitate care in treating the experimental data, which will be 
demonstrated shortly. The important point to be emphasized is that the 
relations given above in Eqs. (7) and (8) hold even in the presence of 
the various mechanisms listed earlier which may occur in modern silicon 
devices', rncluding the effects of carrier drift, position-dependent 
I Lfetimc, and band-edge distortion. 

From Eqs. (5) through (8) , we find that the saturatjon current 
is given by 
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(9) 


^QNO ''' 

and the relaxation time of the diode, including both the emitter 
and base regions, is given by [17] 


T 

II: . 


(Q +0 ) /I 

'^EO ' QNO 


( 10 ) 


The relaxation time characterizes the exponential response of the excess 

minority- carrier charges q„(t) and q„(t). If a steady forvrard current 

is applied and then removed, the excess minority charges q„(t) and q„(t) 

E B 

will decay or relax toward zero exponentially with time as exp (-t/x^) . 

The system of electrons and holes in the quasi-neutral regions relax toward 

the equilibrium state exponentially with time constant x . 

m 

There are six unknowns in Eqs. (9) and (10) : the relaxation time, 

the saturation current, the two effective lifetimes, and the two pre- 
exponential factors defined in Eq. (8) . By methods to be described 
below, the relaxation time and the saturation current can be determined 
from experimental data. If values can be fixed for any two of the 
remaining unknowns, then all six unknowns are determined. In general 
this may require experiments [19] in addition to those described in 
this paper, as well as computer solutions [20] of the differential 
equations underlying the diode characteristics. But, for the diode 
structures studied here, the design justifies additional approximations 

that immediately determine Q and Q_„, simplifying the problem greatly- 

EO BO 

To determine Q^q/ we note that the junctions are deep. Thus the 
carriers injected into the quasi-neutral emitter never reach the highly- 
doped part of the emitter, and the effects of band— edge distortion can 
be neglected. Hence, using this assumption together with the assumption 
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of an exponential doping profile in the emitter , one can show that [21] 




This implies that the electric field due to the gradient of impurity 

concentration, in effect, packs the carriers into a thin region near 

4 

the junction space-charge region. For example, if ~ ' 

this thin region occupies about one-tenth of the thickness of the 
quasi-neutral emitter. Here denotes the doping concentration at 
the emitter edge of the junction space-charge region and is the 

majority carrier concentration at the emitter surface. To a first 
approximation, which neglects the de-ionization of impurity atoms in 
heavily-doped material, equals approximately the surface impurity 

concentration. A suitable value of N to be substituted into Eq. (11) 

A 

can be found from a linearly-graded model of the junction [22] . 

To determine for the diodes under study, we use the observa- 
tion noted in Section III that the densities of the deep energy levels 
in the base show no spatial dependence. This permits the use of the 
widely-known approximation. 


(Aqn,/N^^)/D^T, 


in which is the impurity concentration of the homogeneous base. 

The value of can be found from measuring the dependence of small- 
signal capacitance on reverse bias, a method that is commonly used. 

Hence, for the diode structures designed for this study, determina- 
tion of the contribution of the emitter to the diode characteristics, 
and of the emitter and base lifetimes, requires only that Iqjjq and 
be experimentally determined and then combined with Eqs. (9) through (12) . 
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6.2 Determination of the Saturation Current I 


QNO 

The parameter in model above corresponds to the diffusion 

component (component a) of the diode current described in Section IV. 
Summing this component with the component representing the bulk and 
surface space-charge layers, as in Eq. (1) , yields the total diode 
current. The first step in determining therefore,' is to separate 

the measured static I-V dependence into the two components of Eq. (1) . 
This involves siibtracting the bulk and surface component, which has a 
reciprocal slope ®jj- ^ measured diode current to reveal -the 

diffusion component, which has a reciprocal slope m=l. The procedure 
is shown in Pig. 2. In this figure, the broken line is the extrapolated 
bulk and surface component. The solid curve is the measured I-V 
dependence. Subtraction yields the line that passes "through the dots, 
which is the ideal diffusion component of the diode current and has a 
reciprocal slope m=l. The ideal diffusion component can dominate in 
determining the measured current over a considerable range of current. 
For the example of Fig. 2, the range is about two decades of current. 

At higher currents, the measured current may drop below the ideal com- 
ponent, with m > 1, either because of series resistance or high injec- 
tion or both. 

By the procedure shown in Fig. 2 the saturation current I is 

QNO 

thus simply and unambiguously given by the intersection of the ideal 

diffusion component with the V=0 axis. For the diodes studied here, 

the results thus derived for I_„_ are shown in the first column of 

QNO 

Table III. 


An alternate method exists to determine I 

QNO 

taking the I-V measurement at higher temperatures. 


This simply involves 
which emphasizes the 



Current (A) 



I 1 / ,.,i ! ; L_ 1- —0 

0 0.1 0.2 0.3 0.4 0.5 0.6 


Voltage (V) 

Fig. 2 - Separation of the I-V curve into components. The solid curve 
is the measured I-V dependence. The broken line is extra- 
polated bulk and surface component. Ideal diffusion component 
passes through the dots. 
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dominance of the diffusion component. At a high enough temperature 


the reciprocal slope becomes essentially unity over several decades 
of current, enabling direct extrapolation of the measured current to 
V=0 without the necessity of subtracting and decomposing into com- 
ponents. For diodes having reciprocal slopes within about ten per- 
cent of unity at room temperature, we found that increasing tempera- 
ture by 30°C sufficed to produce unity slope. 

An additional test can be made to determine the voltage range 
over which the diffusion component dominates and thus to lend further 
confidence to the procedure by which is determined. This test 

involves measuring the dependence of forward current on temperature, 
which gives the activation energy. If the diffusion component dominates 
at the forward voltage under study, theory [8] predicts that the activation 
energy will be approximately the band gap. If not, it will be about 
one-half of the band gap. This measurement was made for several diodes 


at voltages in the middle of the voltage range for which the diffusion 
component appeared to dominate. The result deduced from this measurement 
was about 1.1 eV, which is in close agreement with the value of the 
silicon energy gap. 


6.3 Determination of the Relaxation Time 

In general, the relaxation time, x , can be found by measuring 

m 

the diode response to some form of time-varying excitation. Of the 
several methods available to determine the relaxation time this way, 
we used two: the jxinction— current- recovery (JCR) method [15] cind the 

open-circuit- voltage-decay (OCVD) method [16] . 

As the preceding discussion in this section has emphasized, what- 
ever method is used we must confine the measurement to that range of 
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voltages fcr which the diffu-sion component of the diode current 
dominates over the component coming from recombination in the siirface 
,and bulk space-charge layers. In this regard, OCVD has the advantage. 
Fig. 3 illustrates a typical OCVD response. After an initial jump due 
to the cessation of current through the series resistance, the response 
shows a linear variation in time, having a slope, dV/dt. 

The voltage range of this constant slope corresponds to the range 


seen on the static I-V characteristics for which the diffusion component 


dominates. Thus if we assume 




% 

dt 

dt 

^ dt 


in Eq, (6) , then 


kT 1 

^m q dV/dt 


(13) 


which is the relaxation time constant describing the quasi— neutral 
regions. As the voltage decreases beyond this range. Fig. 3 illustrates 
that the slope lessens, a result of the increasing role of the recombina- 
tion component from the bulk and surface space-charge layers. But, by 
focusing on the initial portion of the response and using Eq. (13) , one 
can find the relaxation time which is dominated by the charges stored 
in the quasi-neutral emitter and base regions. This is the relaxation 
time of the model described in Section 6.1. 

In contrast, determination of t by JCR involves some average over 

m 

a rcinge of forward voltage that includes a low-current range where the 

diffusion component is insignificant while the . components from the bulk and 

% 

surface space-charge-layers dominate. 

We used both OCVD and JCR to determine t^, keeping the maximum 
forward voltage applied to the diode below that which would cause high 
injection. Results obtained by both methods are compared in Table II. 
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Time (jjsec) 

Fig. 3 - Typical OCVD response showing the distinctive linear decay in the diffusion range. 



In all discussions to follow in this paper, the values used will be 
those found from the OCVD method, which are the more appropriate values 
for the reasons just stated. 

As was noted in Sec. 6.2, raising the temperature increases the 
dominance of the diffusion component. Thus, in the OCVD response, an 
increase in temperature extends the voltage range for which the initial 
constant slope prevails. This was done for several diodes, and the 
values of the lifetimes were calculated. These lifetimes 
agreed with those determined at room temperature. The results of this 
procedure justifies the neglect of the bracketed components in Eq. (6). 

A variation of the OCVD method could be used which restricts the 
response itself to the range for which the diffusion component is 
dominant. In this variation, the diode is biased in this range and then 
a small pulse of current is superimposed. The voltage response is linear 
with time, showing a single slope. 

Apart from the JCR and OCVD methods, t can also be obtained from 

ra 

measuring immittance of the diode versus frequency [20,22,23]. 

6.4 Emitter and Base Lifetimes and Other Material Parameters 

Combining Iqjjq as determined by the experiments described 

in the preceding sections, and using the theoretical relations described 
in these sections, we obtain the effective emitter and base lifetimes, 

Tg and Tg- In Tables II and III these are listed for each of the four 
substrate doping concentrations being studied. These tables contain 
also the pre- exponential factors Q and Q in Eq. (8) which give the 
charges of the excess minority carriers, Q„ in the emitter and in 
the base. To emphasize the contribution of the emitter to the diode 
characteristics, various ratios involving these material parameters are 
displayed in Tables II and III. 
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TABLE III 


Summary of Material Parameters at 297. 5°K 





VII. DISCUSSION 


The emitter has more influence on the characteristics of the diodes 

having the higher substrate impurity concentrations, which is to be 

expected on qualitative grounds. With regard to the static characteristics, 

Colum G of Table III exhibits how the emitter contribution to the total 

current depends on the substrate concentration. This emitter contribution 

reaches about 40% of the total current for device 2-5-lp5 
1 3 3 

= 6.7x10 cm ) . As column E of Table III shows, this corresponds 
to an emitter current that is about 65% of the base current- This 
emitter contribution for device 2-5-lp5 is more accurately determined 
than the contributions for the diodes having lower substrate concentrations. 
For these devices, the figures listed for the emitter contributions are 
less accurately determined because their calculation involves subtracting 
two nearly equal mambers. 

With regards to the transient characteristics. Column G of Table II 
demonstrates that the base lifetime does not necessarily equal the time 
constant measured by open-circuit-voltage decay. For device 2-5- lp5, 
an error of about forty percent in the calculated base lifetime is intro- 
duced if the emitter contribution to the transient response is neglected. 
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lifetime is considerably lower than the base lifetime in the device 
(2-5-lp5) with the highest substrate doping and the shallowest junc- 
tion. For this device, recombination within the emitter occurs in a 
region of relatively high doping that is relatively close to the 
degenerately-doped surface. 

The results discussed above apply to the particular diodes 
fabricated for this study. They demonstrate quantitatively the im- 
portance of recombination in the emitter to both the static and 
transient characteristics. This demonstration was one purpose of 
this paper. But the main emphasis of the paper is on the method 
which yields the results. For the first time, a method is described 
which separates as well as determines the emitter and base lifetimes 
in a p-n- junction device after the junction has been fabricated. 

The method involves the following measurements: (i) the 

static current- volt age dependence, (ii) the transient open-circuit- 
voltage decay, (iii) the capacitance versus reverse voltage, 'and 
(iv) the thermally-stimulated, voltage- stimulated and light-stimulated 
capacitances [4,5]. These methods of measurement are not new. The 
novelty presented here lies in combining and interpreting the results 
of these various measurements with simple models that are consistent 
with the physics underlying the electrical characteristics of semi- 
conductor p-n- junction diodes. A major result is that the method 
presented here discloses the base and emitter lifetimes concealed in the 
experimental data. 

The design of the diodes used here has permitted the study of the 
emitter and base lifetimes unobscured by the possible effects of energy 
band-edge distortion (gap shrinkage) which may be present in the emitter. 
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The diodes were designed to have deep junctions with moderate impurity 

gradients on siliocn siibstrates having impurity concentrations not 

exceeding 10 atoms/cm . This allows us to calculate the equilibrium 

emitter charge Q defined in Eg. (8) under the assumption of no energy 
EO 

gap shrinkage in the thin portion of the quasi-neutral emitter where many 
minority electrons exist. 

The self consistency of this assumption can be easily demonstrated 
for the diodes studied. From the assumption of negligible gap shrinkage 
it follows that the excess carriers within the emitter are packed by 
the built-in. electric field into a thin region adjoining the junction 
space-charge region. In Section 6.1, we noted that this active portion 
of the quasi-neutral emitter extends from approximately 


" = - 


- — n:; T- ■ - } to X = W„. Hence, the assumption 


MAX' AA E 


of negligible gap shrinkage is self-consistent if N^(x^) .is sufficiently 

small. The assimiption is most questionable for the diode (2-5-lp5) 

having the highest base doping concentration and the shallowest junction. 

For this device, using a gaussian approximation of the doping profile, 

17 3- 

we calculate that N (x ) is about 2x10 atoms/cm . At this impurity 

AA J. , 

concentration, little gap shrinkage can occur [25] ^ which demonstrates 
the self consistency of the assumption. 

But this assumption is questionable for p-n-junction devices 
having shallower junctions, and higher impurity gradients and base doping 
concentrations. For such devices, besides the emitter and base life- 
times, the energy gap shrinkage in the emitter needs to be determined.. 
Variations ■ of the method described here, involving additional measurements, 
have been used to determine the energy-gap shrinkage for p-n-junction 
solar cells and for p-n-junction transistors [26]. In Chapter. 4, we 
describe the application of these methods to n-on-p solar cells. 


62 



appendix a (of CHAPTER II) 


Electron thermal-emission rates from transient- capacitance measure- 
ments for four imperfection levels in the upper half of the bandgap are: 

e^ = (7.7±2.1)10^(T/300)^exp(-93±2)/kgT sec"^ 

e^ = (1.51±l.l)10^^(T/300)^exp(-270±7.4)A„T sec”^ 
n B 

e^ = (8.72±4.5)10^°(T/300)^exp(-285±6.4)/k„T sec“^ 
n B 

= (2.89±1.3)10^^(T/300)^exp{-542±8.6)/k T sec“^ 
n B 

The hole thermal- emission rate for the midgap level, obtained from 
measurement of diode reverse leakage current [4] , is given by 

= (1.2±0.46)10^ (T/300)\xp(-593±9.5) AgT sec”^. 
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CHaPTER IV. 

METHODS FOR DETERMINING THE EFFECTIVE GAP SHRINKAGE 
(MINORITY-CARRIER STORAGE) AND EFFECTIVE LIFETIME IN THE EMITTER 

I. INTRODUCTION 

The electrical characteristics of p-n-jiinction solar cells depend 
upon the recombination, generation, and transport properties of the mobile 
electrons and holes within the device material. These properties were first • , 

described by Shockley [1] in his ideal diode theory, which treated the 
relevant material parameters as constants. But within the emitter of a • ' 

modern silicon cell the behavior of these properties can be more complex 
than that assumed in the ideal diode theory because of the presence of band- ! 
edge distortion (energy-gap shrinkage) and of position-dependent recombina- • 
tion rates [2] . As was discussed in Chapters I and II, these two mechanisms < 
have been suggested as probable reasons for certain discrepancies between * 

the theoretical and experimental performance, particularly the discrepancy of 
loo millivolts in the open-circuit voltage seen in diffused-emitter silicon 
solar cells made with low-resistivity substrates [3] . Previously the physical 
mechanisms occurring in the emitter have been the subject of speculation and 
theoretical studies; no means has been available to investigate and character- 
ize these mechanisms by experiment. 

This chapter describes an experimentally-based methodology for determining 
the effective lifetime and the effective gap shrinkage (minority-carrier storage) 
in the emitter of p-n- junction solar cells. The methodology pertains to a solar- 
cell structure after the junction has been fabricated, its use distinguishes 
the emitter region from the base region, yielding material parameters that apply 
•to each of these two regions. Thus it enables a determination of the relative 
roles of the emitter and base in governing solar-cell performance, and an assess- 
ment of the importance of gap shrinkage relative to that of large recombination 
rates in the highly-doped emitter. 
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This chapter will emphasize the methodology — that is, the strategy of 
the procedures used and the physical principles underlying them. In the 
organization of the chapter, some of the theoretical grounds for the method- 
ology are first described. This leads naturally then to the experimental 
measurements used and the physical interpretation of the data. To illustrate 
the procedures, they are applied to an n'*'-p- junction solar cell having a 
siibstrate resistivity of about 0.1 ohm-cm. 

II. A MODEL RELATING STATIC 2^D DYNAMIC RESPONSE 

To provide a basis for studying the material parameters of the emitter 
and the base, a model of diode behavior is needed that is general enough 
to apply: 

(a) in the presence of high, position-dependent recombination 
rates and of gap shrinkage; 

(b) for any doping profile associated with the p-n-junction. 

Recently such a model was proposed [4] . This model was outlined in Chapter III 
but we sketch its development here again to emphasize certain of its features 
that bear on the methodology and to enable this chapter to be read indepen- 
dently of Chapter III. 

As is discussed in Chapter I and in the Appendix, an understanding of ' 
cell behavior in the dark is of central importance to the discrepancy 
between theoretical and experimental values of V seen in low-resistivity 
silicon cells. Thus, consider a p-n-junction diode in the dark. In the 
ideal model of Shockley, the diode is divided into two types or regions: 

(a) the quasi-neutral emitter and base; and 

(b) the junction space-charge or transition region. 

In a practical silicon diode, other space-charge regions may also exist, 
at the oxide-silicon interface and at the current contacts [5] . 
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We now write the static (or d.c.) current, I', as the components 


from the quasi-neutral emitter and base regions added to the components 

coming from the bulk and surface space-charge regions. By the princi- - 
pies of charge control [6] , 


Here Q and Q designate the magnitude of the charge of the excess 

minority carriers in the quasi-neutral emitter and base, and SO denotes 

^SC 

the effective moble charges stored in the bulk and the surface space- 

charge regions. The parameters x_, x„, and x__ are the relevant charge- 

£ B SC 

control time constants or relaxation times. 

To describe the current i ( t) under dynamic conditions , we write 


i(t) = ^ + ^ 

■^E "^B 


dq 


^sc 

'^sc 


dq 


+ 


SC 


dt 


( 2 ) 


which follows from the quasi-static approximation used in charge-control 
analysis [6] . In Eqs. (1) and (2) , the components of current from the 
space-charge regions are bracketed for emphasis. 

The methodology presented in this paper will emphasize the material 
parameters describing the quasi- neutral emitter and base. Hence the 
components of the current arising from the space-charge regions will be 
neglected in Eqs. (1) and (2) . This approximation simplifies the model. 
But it requires that care be taken in experimental measurements; either 
conditions must be designed to make these space- charge-region components 
negligible or they must be removed from the experimental data to reveal 
the behavior coming from the quasi-neutral regions. In Sections III and 
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IV of this paper, we will describe specific procedures for dealing with 
the space-charge-region components and for dealing with other assumptions 
to bo made in developing the model for diode behavior. 

[ r this approximation is combined with the assumptions of low in- 
jection and negligible series resistance, the currents and the excess 
charges can be written simply: 


i = IgjjQ[exp(qvAT) - 1] (3) 


Here v = v(t) .designates the total time- dependent voltage at the terminals 
of the diode. The subscript QN emphasizes that the saturation current 
^QNO from processes occurring within the quasi-neutral regions. 

The additional assximptions needed to yield Eqs. (3) and (4) also require 
special care, to be discussed, in the design of experiments and in the 
interpretation of the data. 

Combining Eqs. (1) through (4) gives 

(5) 


= O /t + Q /t 
QNO EO B BO B 


,\ 


for the saturation current. By setting i(t) =0 in (2) and solving the 
resultant differential equation, we find [4] 


^QN ^^EO ^ 


( 6 ) 


for the relaxation time of the model. This relaxation time determines 
the natural (force- free) behavior after the device is perburbed by an 
applied force, such as by an impulse current or a voltage step. It is 
an important parameter in determining the response of the diode to any 

* t 

time-varying (current) driving force. 

Although simple, Eqs. (5) and (6), prove to be powerful relations 
for determining parameters of the quasi- neutral material. As will be 
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seen, the charge is closely related to gap shrinkage occurring in 
the emitter; and is the effective lifetime of the minority carriers 
in the emitter, which relates to the position-dependent recombination 
rates that may be present there. Hence, it becomes of interest to see 
how and — and their counterparts for the base and x^ — can be 

determined from experimental data. 

The possibility of doing ‘this depends on which of the six unknowns 
appearing in Eqs. (5) and (6) can be measured and which can be calculated 
by auxiliary relations. There are many possible measurements that apply 
to this problem, and some of these will be discussed in Section III. 

The measurements discussed in Section III will pertain to the determine- 

^bn' ®bo- ''b- 

Theory offers, in addition, certain auxiliary relations. First, if 
the minority-carrier lifetime and the net doping concentration both stay 
constant with position in the base, then from the conventional theory for 
low injection, we have for the base region. 


«B0 ' ® ■ 


(7) 


where L = /d x is the minority-carrier diffusion length in the base, n. 
is the intrinsic carrier concentration, and N is the doping impurity 
concentration in the base for the n^-p structure treated here. The spa- 
tial constancy of and x^ can be assessed by the capacitance and 
transient-capacitance measurements indicated in Section III. 

Besides Eg. (7), there is a second ' theoretical relation of interest. 
Let designate the charge in the absence of gap shrinkage in the 
emitter. By a simple extension of conventional theory, a relation for 

can be found [7] . Pig. 1 -illustrates the physics underlying this 
EO 
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relation. The gradient of the impurity concentration in the emitter 
produces a built-in electric field that packs the injected minority 
carriers into a thin region near the edge of the junction space-charge 
region: Integrating the density P (x) of the minority holes yields [7] 



qn./N(W ) 

3 . 


] 


f % ^ 


( 8 ) 


in which is the maximum majority-carrier concentration and N(W^) is 

the majority-carrier concentration at the emitter edge of the junction 
space- charge region. From Eg. (8) , to a first approximation, the minor- 
ity carriers can be viewed as restricted to the region,- W 

s 

which can be regarded as the active portio n of 

the emitter of the dark diode in the absence of gap shrinkage. This point 

of view is emphasized by the dashed rectangle drawn in Fig. 1. 

In eq. (8) , x = W corresponds to the emitter edge of the junction 

space-charge region. If the doping profile is assumed to be distributed 

spatially as any monotonically decreasing function, such as a Gaussian 

or an exponential function, then will also be the thickness of the 

quasi-neutral emitter and will’ be the majority- carrier concentration 

MAX 

at the surface. But if the profile is practically flat over' part of the quasi- 
neutral emitter, then will stand for the distance from the edge of the junction 
space-charge region to -the plane where the doping concentration starts to 

rapidly decline. To get an accurate estimate of an accurate estimate 

EO 

of the built-in electric field near the junction space-charge region is 

required. The meaning of W just discussed helps provide this estimate. 

E 

It applies unless the doping concentration is practically flat over the 
entire quasi-neutral emitter, in which case the built-in elec-tric field 
i ‘5 nearly zero and Eq. (8) is invalid. 
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Pig. 1 Tile distribution P(x) of the excess minority carriers (holes) in 
the quasi-neutral emitter of a dark solar cell. The surface is 
at X = 0 and the edge of the junction space-charge region is at 
X = Wg. The area of the dashed rectangle equals the area under 
.the curve/ P(x). 
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Ecjuation (8), defining Q°_, has several uses. First, it pe 2 nnits an 

HU 

approximate calculation of the emitter doping density at 

, which is the maximum density in the active portion 

18 3 

of the emitter. If this density falls below about 10 atoms/cm , the gap 
shrinkage in the active portion is negligible [8] . Hence, as a trial 
calculation, the gap shrinkage can be assumed to be negligible and 
can be calculated from Eq. (8) . This assumption can then iiranediately be 
tested for self-consistency , in the manner just described. For devices 
having junctions deep enough as well as small enough impurity gradients 
and base doping concentrations, this assumption will be self-consistent. 

Thus calculated from Eq. (8) can be combined with Eqs. (5) and (6) and 
with the results of various measurements to yield values of t and t , the 

H 8 

emitter and base effective lifetimes. This has been done as part of a 
study on a set of specially designed diodes [9] • 

A second use of Eq. (8) is as a defining relation for an effective 
gap shrinkage Ae^ in the emitter. The definition follows directly from the 
meaning of gap shrinkage employed in [2] and discussed in Chapter II; 

' exp(AE^AT) (9) 

This definition of Ae agrees exactly with the meaning of Ae„(x) discussed 

G 

in [2] for the case that the net doping concentration in the quasi-neutral emit- 
ter is flat and the emitter minority- carrier lifetime is independent of 
Ae^(x). For other cases it gives a weighted average of the energy-gap 
shrinkage over the thickness of the emitter. 


JL ^ 

/M(W„) 
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III. MEASUREMENTS AND INTERPRETATION 


This section treats various measurements -that lead to values for 

'^T,’ which axe parameters appearing in Eqs. (5) and 

o du QN 

(6) of the model described in Section II. That model assumes the 
dominance of processes occiirring in the quasi-neutral emitter and base 
regions over those located in the junction space-charge region or at 
the surface. Hence the discussion here will deal with assuring this 
dominance or, alternatively, with removing from the experimental data 
that part which belongs to processes in the junction space-charge region 
and at the surface. A key to appropriate design of eiqjeriments and proper 
treatment of data is the careful linking of that data to the physics 
governing the observed behavior. 

3.1 Device Structures 

The discussion of this paper will limit consideration to high- 
efficiency junction soiar cells for which processes occurring in the 
quasi-neutral emitter contribute significantly to the cell performance. 

Hence the device structures of interest are made of single-crystal semi- 
conductors, such as silicon. For certain of the measurements - for example, 
open-circuit voltage and short-circuit current ~ an entire solar cell 
may be used. For transient or small-signal measurements, a smaller area 
is desirable, and the structure used may be part of an entire solar cell 
or it may be a junction diode specially designed for the purpose. In the 
measurement of the properties of a dark cell, a diode whose metallization 

geometry contains an MOS guard ring offers advantages, particularly for 
4 - 

n -p diodes. A voltage applied to the guard-ring can effectively eliminate 
the component of current arising from recombination at the surface. 
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3.2 Base T- and Q„ 

In the special case of spatially- independent and N^,-Eq. (7) 

links Tg and together with the diffusion length = v^D^x^ of the 
minority electrons. Both and x^ may be explored by various 
measTirements of the reverse-biased capacitance C- For example, 
measuring 1/C versus V, the reverse bias, determines the value of 
N' if it is independent of position. By plotting the same data as 

3 

l/C versus V for small V, one can infer the impurity grade constant 
^ of the junction, which enters in the determination of N(W^) 
appearing in Eq. (8) . 

By various transient-capacitance methods [10] the value of x^ atd 
its spatial dependence can be found. These methods include: 

(a) Thermally-Stimulated Capacitance (TSCAP) , which 
involves measuring high-frequency capacitance versus, 
temperature; 

(b) Voltage-Stimulated Capacitance (VSCAP) , which involves 
measuring high-frequency capacitance versus time after 
a bias voltage change; and 

(c) Light-Stimulated Capacitance (LSCAP) , which involves 
measuring high-frequency capacitance versus time after 
a change in the illumination. 

In many junction solar cells x^ and will be independent of 

X, validating Eq. (7). For such cells, the diffusion length L , and 
• n 

hence the lifetime x^, can be independently determined by a technique 
involving the use of X-rays. In this technique [11] the short-circuit 
current generated by a uniformly- absorbed beam of 250 KeV X-rays is 
collected by the junction. This current is nearly proportional to L^, 
and a calibrated relation exists [111 for determining L . 
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3.3 Saturation Current ^ 

QNO 

The parameter ^qj^q can be inferred either from the static I-V 
characteristic of a dark diode or from the static (photo- current 

versus photo-voltage) response of the diode to various levels of illumina- 
tion . The key to inferring from the data lies in the proper 

linking of that data with the basic mechanisms producing it. 

3 . 3 (A) Dark current measurement ; The data are taken and plotted as 
In I versus V. The tmderlying physics requires that the data follow 
the relation, 


I = igjjjQ[exp(qV/kT)-l] + I^Q[exp{qV/m^kT)-l] , 


( 10 ) 


in which the first term describes the component of the current that 
comes from the quasi-neutral regions and the second term describes the 
component coming from the bulk [12] and surface space- charge regions [5] . 

As a first step toward determining observed dependence 

of J-n I versus qV/kT is resolved into the two components of Eq. (10) . 

This can be done graphically by subtracting the component involving 
reciprocal slope m^, which dominates at low current, from the observed ‘ 
dependence to reveal the component with unity reciprocal slope. This 
component intersects the ^^n I axis at An Iqjjq' according to Eq. (10), and 
the saturation current is thus determined. 

The physics underlying Eq. (10) does not include the effects of series 
resistance or high injection, both of which can bend the curve of An I 
versus V at high currents, increasing the reciprocal slope. For many highly-doped 
cells this effect does not obscure the two components of Eq. (10) because 
the bending occurs at currents several decades higher than 
the low-current range for which Eq. (10) applies. If a clear separation 


76 



does not exist, Eq. (10) must be modified to include the relevant 
physics [13] . 

3.3(B) Photo-current measurement ; The short-circuit current and 
the open-circuit voltage are measured at various levels of 
illumination, and the curve of £.n versus qV^^/kT is plotted. 

This curve follows the relation, 

^SC "" ~ exp (*5130^3 At)] + I^^[exp(qV^^/m^kT) - exp(qIg^Rg/m^kT)] 

( 11 ) 

which is consistent with the underlying physics. For completeness, 
the effect of a series resistance R is included, but Eq. (11) must be 
modified further if the effects of high injection should prove to be 
significant in a particular cell. As in the treatment of the dark-current 
data, the measured dependence is resolved into the two components given 
in Eq. (11) . Then is the intercept on the £n I axis made by the 

component with vinity reciprocal slope. 

3.3(c) Temperature Dependence ; As an alternate method for determining 

^QNO ^ check against the results of the methods just described, 

the current-voltage dependence can be measured at temperatures higher 

than 300®K. The component of current coming from the quasi-neutral regions 

2 

is approximately proportional to n^, which in turn is proportional to 

exp(-Eg/kT), where E^ is the band gap. In contrast, the components 

arising from processes other than those located in the quasi-neutral regions are 

2 

approximately proportional to n^ rather than n^ and thus to exp(-E^/2kT) 
rather than exp(-Eg/kT). Accordingly, as temperature rises the com- 
ponent of current from the quasi-neutral regions will tend to become 
more dominant and the measured dependence will tend more toward unity 
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initial charges residing within the device, and 
then a reverse current is applied to speed the decay 
of these charges. 

(b) Open-circuit voltage decay [15] , or OCVD, in which 
the initial charges are set as in JCR but then are 
allowed to relax in the open-circuited condition 
with no reverse current applied. 

(c) Admittance vs. frequency [16] , in which a small 
sinusoidal excitation imposed on the bias condition 
is applied and the real and imaginary parts of the 
complex admittance are ni^asured as functions of 
frequency and bias. 

The relaxation time governs the force-free response of the 
model of Section II, in which the charge of the mobile carriers within the 
quasi-neutral regions is assumed dominant over that within the bulk and 
surface space-charge regions. Therefore, if the quasi-neutral regions 
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dominate, the OCVD response directly gives OCVD method has 

an advantage over the JCR method because JCR necessarily involves 
an excursion over the entire range of forward voltages including 
that range in silicon diodes for which the space-charge region will 
dominate in producing the response. Hence the relaxation time in- 
ferred from JCR data will differ from . In contrast, the initial 

QN 

linear portion of the OCVD response, from which the relaxation time 

is deduced, can derive from processes located in the quasi-neutral 

regions, and thus can determine . 

QN 

Whether the processes in the quasi-neutral regions dominate in 
producing the dynamic response is disclosed by the dependence on 
bias voltage seen in the real and the imaginary parts of the measured 
admittance [17] • Increasing the temperature will increase the dominance 
of the quasi-neutral regions , for much the same reasons as 
in Section 3.3(c). 

3.5 Other Measurements 

The measurements discussed above are electrical. Other measurements, 
electrical and non-electrical, could be taken with possible advantages. The 
of measurements given here is meant to be suggestive, not exhaustive. 

3.6 Internal Consistency and Measurement Sensitivity 

The experimental results yielded by these various measurements 
overlap. Because they are not mutually exclusive, some checks can be 
made for internal consistency. For example, from Eqs. (5) and (6) 
alone it follows that: 


4 ^ 5 /Q = - 1 

^EO '^BO Q ^ 

^BO 


( 12 ) 


TbAe = 


(13) 


79 



where A = Q /O = [exp(AE„/kT) is a measure of the gap shrinkage. 

EO “BO G EO BO 

Alternatively, expressions for A and t„/t_ can be developed that 

S £| 

relate to V and I the measured open-circuit voltage and short- 
circuit current of the illiminated cell. To develop these, we use [18] 


“ W^SC • 


(14) 


Here is the open-circuit voltage that would occur in a solar cell 

if the dark current resulted from recombination only in the quasi-neutral 
regions. In other words, is the open-circuit voltage of an 

idealized solar cell for which the second term in Eq. (10) , representing 
space-charge-layer recombination, is negligible when I = To correct 

for the presence of recombination in the bulk and surface space-charge 
layers, the expression 


/kT = qV„„AT + in 
OCqn OC 


SC 




(15) 


relates the measured voltage V„_ to the idealized voltage . This 

expression follows directly from inspection of the graph of the two 
components of I in Eq. (10) . The two components are determined from the 
measured cur rent- voltage dependence as explained in Section 3.3(A). 


The voltage accounts for recombination in both the quasi- 

neutral base and emittcar regions. We identify 


^Cqn 


kT 


r I 


in 


SC 




(16) 


as the open-circuit voltage if recombination in the base alone were 
important, and 
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(17) 


6V 


oc 


= V® 
OCqn 


- V 


OCqn 


as the shift in open-circuit voltage resulting from significant recom- 
bination in the emitter. Note that is positive. Combining Eqs. (5) 

and (6) with (14) through (17) yields 


exp (qSV /kT) 

A = -22 (18) 

Tg/Xg = (1/A) [exp(q5V^^/kT) - 1] , (19) 

which are the desired expressions. 

Eqs. (12) and (18) provide two alternative calculations of A, and 
hence of gap shrinkage. Eqs. (13) and (19) provide two alternative 

means for calculating the emitter lifetime. None of these four 
relations depends strongly on the intrinsic concentration, n^. This is 
desirable because n^ and its temperature dependence is not accurately 
known at present. In contrast to Eqs. (18) and (19) , note that the use of 
Eqs. (12) and (13) does not necessarily require measurements \mder illumina- 
tion. This indicates that the two methods of calculation are, to 
some -degree, independent. As part of the methodology we require for any 
device under study that the calculations of A by Eqs. (12) and (18) agree 
satisfactorily, and that the calculations of x by Eqs. (13) and (19) like- 

Cl 

wise agree. Expressed in general terms, internal consistency is required. 

Eqs. (12) and (18) in effect express the charge as the difference 
between the charge and the total (excess minority-carrier) charge, 

For some devices — for example, those in which gap shrinkage is small — these 
two equations call for the determination of Q by subtracting two numbers 

£iU 

that nearly are equal. Because some imcertainty will always exist about 
the meastired values used in these calculations, the method will give no 
accuracy in the calculation of A), for these cases. 
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The measurement sensitivity of the effective gap shrinkage and 


lifetime of the emitter thus depends on the accuracy with which each 

of the measured quantities— and are 

determined. Let A . be the minimum value of A that can be determined 
min 

for a certain precision of the measurements. (For example, if = 1/ 

it is assumed that Q can be determined from Eqs. (12) or (18) only 

EO 

if Q >0 .) Then, from Eq. (8) and the definition of A, the minimum 

EO - BO 

determinable value of the gap shrinkage Ae . is 

G « 

min 


Ae„ = kT£n[A . 

G , min^BO ^EO 

min 


( 20 ) 


= k«n{ «.nW^(Wj,)) 

This is defined as the measurement sensitivity for the gap shrinkage 
of a given device. Similarly, combining Eq. (20) with Eqs. (13) or 
(19) , we can define the corresponding measurement sensitivity for the 
emitter lifetime . 

3.7 Path Toward Greater Measurement Sensitivity 

From Eg. (20) we note that increasing A . L /W_ will give the greatest 
^ min n E 

increase in measurement sensitivity. Eqs. (12) ^d (18) imply that the quantJ 

A . is set by the precision with which T , L , N , V , and I 

min ^ Qjj QjjQ n AA OC sc. 

are measured. The emitter thickness W can be purposefully increased 

£j 

in cells especially designed to explore the gap shrinkage experimentally. 

Some evidence exists that the diffusion length decreases rapidly with 

18 —3 

increasing doping concentrations in the base for >10 cm 
Hence, increasing the conductivity of the base region could much improve 
the sensitivity. Alternatively, the diffusion length can in effect be shortei 
by structural and electrical means through adding a collector region, making 
a transistor- like structure. The third terminal that can be contacted 
to such structures also provides for a direct and accurate measurement of 
Tg, from which and the gap shrinkage can be deduced [19] . 
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IV. ILLUSTRATIVE EXAMPLE 


4. 1 Device 'Fabrication 

Six n -p silicon solar cells shown in Fig. 2 were constructed at 

NASA Lewis Research Center using a standard processing technique- The 

1x2 era, vacuum- float- zone, <111> orientation, 0-1 ohm-cin {approximately) , 

p-type silicon wafers were first chemically polished in an HF-HNO^ - glacial 

acetic acid solution. The wafers were then phosphorus diffused in a 

POCl^ atmosphere for 45 minutes at a temperature of 880"C to form the 

n -p junction. The junction depth is about 0.25ym and the 

20 -3 

surface concentration is about 2x10 cm . After removal of the 

phosphorus oxide glass with HF, silver- aluminum contacts were thermally 

evaporated to form the large area contact to the p-type silicon substrate 

and the gridded contact to the thin diffused layer. Contact adherence 

was insured by sintering the cells in argon at 550°C for several minutes. 

For the measurements involving illumination, the entire solar-cell 

-4 2 

structure was used. Mesa structures of small area (3.5x10 cm ) were 
etched for use in the dark-current and the OCVD measurements. 

4.2 Measurements and Calculations 

To illustrate the methodology, we describe here the measurements and cal 

tions made for a typical sample. The data were taken at temperature T = 29,5. 

First the material parameters for the base are found. From the 

X-ray method L^ = 82 microns, and from measuring the reverse-biased 

17 -3 

capacitance the doping concentration is deduced to be N = 2x10 cm 

2 

Prom standard references [20] , this implies a diffusivity, = 10 cm /sec, 

-19 

and hence a lifetime, t = 6.7 microseconds. From Eq. (7) , Q = 5.4x10 C/< 

BO 

o —22 2 

At this point we also calculate == 2.4x10 C/cm , from Eq. (8). 
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The remaining parameters are measured as described in Section III. 

By the OCVD method, = 0.67 usee. The current-voltage measurements made 

-12 2 

on a dark diode give = 1.1x10 A/cm . This value agrees well 

with that determined from illuminated current-voltage characteristics 

measured using a calculator-controlled-data-acquisition system. A 

xenon- light solar simulator used for the AMO {outer-space condition) 

2 

measurements yields V = 0.600 volts and I = 23 mA/cm . The data 

also cfive 6V = 67mV as the voltage shift defined in Eq. (17) . 

OC 

From these data, the calculations of the gap shrinkage and the 

emitter lifetime can be made. First, using data secured only from the 

response of the dark device, we find from Eq. (12) that A = 0.38 

-19 2 

(implying = 2x10 C/cra and AE^ = 170meV) and from Eq. (13) that 
T /t = 34 (implying T = 195nsec) . Second, using data secured by 
measurements using illumination, we find from Eqs. (18) and (19) nearly 

the same values for A, AE^, and t^. 

4.3 Interpretation of Results 

The results of Eqs. (12) and (18) for the effective gap shrinkage 
agree well, so do the results of Eqs. (13) and (19) for the effective emitter 
lifetime. Moreover, alternate methods used for determining 'fgN 

show good agreement. Hence, as discussed in Section 3.6, the requirement 
of internal consistency among the measurements is met. 


Consider now the measurement sensitivity, treated in Section 3.6. 


If 5% accuracy is assumed in the measurement of each parameter 

and Qgg) in Eq. (12) , then the minimum detectable value of A is 


A . - 0.15. From Eq. (20), this corresponds to Ae_ 

min G . 


147meV 


min 

the minumum detectable value of the gap shrinkage for the particular 

solar-cell structures treated. The results obtained here are A = 0.38 

and AE = 170meV which lie within this assumed measurement sensitivity. 
G 
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Greater measurement sensitivity is desirable, however, and to increase the 
sensitivity we have undertaken all of the three approaches indicated in 
Section 3.7. The results of these studies and a detailed consideration of 
the measurement accuracy and sensitivity will be reported elsewhere. 

We emphasize that the methodology described here determines and 

the effective lifetime, and the minority-carrier charge stored in the emit 

These two variables are fundamental in the sense that the values determined 

for them are independent of the specific physical mechanisms assumed in the 

physical model describing the degenerately doped emitter. If one uses a 

physical model for the emitter that assumes the presence of energy gap 

shrinkage, then a simple measure of the effective gas shrinkage can be , 

calculated from Eq. (9). But the fact that the experimentally determined 

value of much exceeds the value Q° calculated assuming the absence 

of gap shrinkage, as is true for the cell studied here, does not itself s 

prove the existence of sizable energy— gap shrinkage. Rather it can only 

show consistency with a physical model containing energy-gap shrinkage. 

A less fundamental problem in interpretation also deserves mention. 

If Eq. (9) is used to relate Q to Ae , then Q° must be calculated. That 

EO , ^ EO 

requires knowledge of the doping concentration profile quasi- 

neutral emitter and, for a shallow- junction cell, it may require knowledge al: 
of the effective surface recombination velocity at the front surface. 

If (x) has a large enough gradient over a considerable part of the 

quasi-neutral emitter, the resulting built-in electric field causes Q° 

EO 

to be practically independent of the lifetime 'the surface 

recombination velocity s. Then Q° is given by Eq- (8) and the value for 

I t i . 

I 

Ae stated above (170 meV) holds. But if (x) were nearly flat over prac- 
o DD if T 

tically the entire quasi-neutral emitter, then the calculation of Q° would 

EO 

require t (x) and s. For this case, if the emitter were transit-time limited 

It 
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rather than lifetime limited, Q° could be considerably larger than the 

EO 

value calculated above and AE could be considerably smaller. 

G 


87 



REFERENCES FOR CHAPTER IV 


1. W. Shockley, "The theory of p-n junctions in semiconductors and p-n 
junction transistors," Bell Sys. Tech. J. , vol. 28, pp. 435-489, July 
1949. 

2. F.A. Lindholm and C.T. Sah, “Fundamental electronic mechanisms limiting 
the performance of solar cells," IEEE Trans. Electron Devices , this 
issue. 

3. H.W. Brandhorst, Jr., "Current status of silicon solar cell technology," 
Technical Digest, 1975 Int. Electron Devices Meeting (75CH1023-1 ED) , 
pp. 331-338, Dec. 1975. 

4. F.A. Lindholm and C.T. Sah, "Normal modes of semiconductor pn junction 
devices for material-parameter determination," J. Appl. Phys., vol. 47, 
pp. 4203-4205, Sept. 1976. 

5. C.T. Sah, "Effect of surface recombination and channel on p-n junction 
and transistor characteristics," IRE Trans, on Electron Devices , vol. 

ED-9, pp. 94-108, Jan. 1962. 

6. E.O. Johnson and A. Rose, "Simple general analysis of amplifier devices 
with emitter, control, and collector functions," Proc. IRE , vol. 47, 
pp. 407-418, March 1959. 

R.D. Middlebrook, "A modern approach to semiconductor and vacuum device 
theory," IEEE Proc. , vol. 106B, suppl. 17, pp. 887-902, March 1960. 

P.E. Gray, D. DeWitt, A.R. Boothroyd, and J.F. Gibbons, Physical Elec- 
tronic and Circuit Models of Transistors , Wiley, New York, 1964. 

7. F.A. Lindholm, "Simple expressions for the minority- carrier charge in 
the emitter of p-n- junction devices," to be published. 

8. F.A. Lindholm, S.S. Li, and C.T. Sah, "Fundamental limitations imposed 
by high doping on the performance of pn junction silicon solar cells," 
Record 11th Photovoltaic Specialists Conf . , (75CH0948-OED) , pp. 3-12, 

1975. 

9. A. Neugroschel, F.A. Lindholm, and C.T. Sah, "A method for determining 
the emitter and base lifetimes in p-n- junction diodes," IEEE Trans . 
Electron Devices , submitted for publication. 

10. C.T. Sah, L. Forbes, L.L. Rosier, and A.F. Tasch, Jr., Solid-State 
Electronics , vol. 13, pp. 759-738, June 1970. See also reviews treating 
these and more recent methods in C.T. Sah, "Bulk and interface imperfec- 
tions in semiconductors," Solid-State Electronics , in press, and 

C.T. Sah, "Detection of recombination centers in solar cells from junction 
capacitance transients, IEEE Trans. Electron Devices , this issue. 

11. J.H. Lamneck, Jr., "Diffusion lengths in silicon obtained by an X-ray 
method," NASA TM X-1894, Oct. 1969. 


88 



12. C.T. Sah, R.N. Noyce, and W. Shockley, "Carrier generation and recombina- 
ation in p-n junctions and p-n junction characteristics," Proc. IRE, vol. 

45, pp. 1228-1243, Sept. 1957. 

13. M. Wolf and H. Rauschenbach, "Series resistance effects on solar cell 
measurements," Advanced Energy Conversion , vol. 3, pp. 455-479, Apr. -June 
1963. 

14. R.H. Kingston, "Switching time in junction diodes and junction transistors," 
Proc. IRE , vol. 42, pp. 829-834, May 1954. 

H.J. Kuno, "Analysis and characterization of p-n junction diode switching," 
IEEE Trans, on Electron Devices , vol. ED-11, pp. 8-14, Jan. 1964. 

15. S.R. Leder handler and L.J. Giacoletto, "Measurement of minority carrier 
lifetime and surface effects in junction devices," Proc . IRE , vol. 43, 
pp. 477-483, April 1955. 

16. W. Shockley, Electrons and Holes in Semiconductors , pp. 313-318, Van. 
Nostrand, Co., New York, 1950. 

H. Maes and C.T. Sah, "Application of the equivalent circuit model for 
semiconductors to the study of Au-doped p-n junctions under forward bias," 
IEEE Trans. Electron Devices , Vol. ED-23, pp. 1131-1143, Oct. 1976. 

17. P.A. Lindholm and M.B. Callaway, to be published; also C.T. Sah, "Effects 
of Electrons and Holes on the Transition Layer Characteristics of Linearly- 
Graded p-n Junctions," Proc . IRE , vol. 49, pp. 603-618, March 1961. 

18. M.B. Prince, "Silicon Solar energy converters," J. Appl. Phys ., vol. 26, 
pp. 534-540, May 1955. 

19. A. Neugroschel, F.A. Lindholm, and C.T. Sah, "Experimental determination 
of emitter lifetime and band-edge distortion in transistors," Electronic 
Letters , to be published. 

20. M.B. Prince, "Drift mobilities in semiconductors. II. Silicon," Phys . 

Rev . , vol. 93, pp. 1204-1206, March 1954. 

S.M. Sze and J.C. Irvin, "Resistivity, mobility, and impurity levels in 
• GaAs, Ge, and Si at 300°K, " Solid-State Electron., vol. 11, pp. 599-609, 

1968. 


89 



CHAPTER V. DISCUSSION 


This research program concentrates on the open-circuit voltage 
in (n-on-p) single-crystal silicon' p-n junction solar cells of 0.1 
ohm-cm substrate resistivity. The purpose is to explain, in funda- 
mental terms,, the reasons for the discrepancy that exists between 
observed values of the open-circuit voltage and those predicted by 
traditional theory — a discrepancy of about 100 mV. 

In Chapter I, we identified anomalous values of the stored minority- 
carrier charge and of the phenomenological lifetime in fhe 
emitter of the solar cell as the probable origins of this discrepancy. 

In Chapter II, proceeding from a theoretical standpoint, we identified 
high recombination rates and effective gap shrinkage in the quasi- 
neutral emitter region as the two main fundamental electronic mecha- 
nisms responsible for the discrepancy in a one-dimensional model of 
the solar cell, and noted the possible role of areal Inhomogeneity 
in contributing to the discrepancy. We examined the many mechanisms 
that can give rise to high recombination rates (low lifetime) and to 
gap shrinkage (large minority-carrier storage) . The lack of relevant 
experimental data and the inadequacies of the existing theoretical 
models prevented our reaching firm conclusions about the role of low 

lifetimes compared with that of gap shrinkage in producing the V 

oc 

discrepancy. Thus, we turned to the problem of experimentally de- 
termining the effective lifetime and effective gap shrinkage in the 
emitter. Our findings are described in Chapters III and IV. 

Although our objective in the research program focusses on n-on-p 
solar cells, the technical findings reported in Chapter III come from 
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measurements made on p -n diodes. We chose a p -n structure for our 
initial experiments for three reasons; 

(1) The charge in the oxide covering the semiconductor surface 
invariably is positive. Hence surface channels and con- 
sequent recombination that obscure the bulk recombination 
being studied tend to occur more in n^-p devices than 

in p^-n devices. This surface channeling can be pinched- 
off by application of guard-ring voltages, but the study of 
p -n devices remains easier than that of n -p devices • 

(2) Schottky-barrier metal- semi conductor junctions must be made 
to find the density of imperfection centers in the substrate 
prior to junction fabrication. These enable study of the 
imperfections caused by the processing used in junction 
fabrication. Schottky barriers can be easily made on n 
material, using aluminum as the metal. For p material, 

the choice of metal becomes a problem of finding one that 
will yield a sizeable barrier height yet remain stable. 
Although not insurmountable, this problem makes the study 
of p^-n cells easier. 

(3) Transient-capacitance and related measurements applied to 
a p^-n cell readily and simply yield information about the 
imperfection centers that lie in the upper half of the 

band gap. These centers control the minority carrier 

* "S' 

lifetxme in the substrate of an n -p cell, and thus study 

of p -n structures yields information of direct relevance 

to the characterization of n-on-p solar cells . 
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of these reasons, the third was the most important in our choice 
to study, n -p devices first- The motivation was to get relevant in- 
formation as quickly and easily as possible while developing under- 
standing and experimental procedures that could be later applied 
(Chapter IV) to n^-p structures - 

Additionally, to get concrete results fast, and being aware of 
the potential problems in sensitivity of measurements done directly 
on material of 0.1 ohm-cm resistivity made at low diffusion tempera- 
tures, we fabricated structures of different doping densities in the 
base (10^^ to 10^^ cm using different diffusant temperatures (900° 
to 1200°C) . Although temperatures above 900®C are rarely used in 
solar-cell fabrication, the high-temperature treatment tends to create 
many imperfections that can be easily measured even on a highly-doped 
substrate. 

In Chapter III, a method was described that provides an experimental 
means for the first time to separate and determine the emitter and base ' 
lifetimes in a p-n junction diode after the junction has been fabri- 
cated. In the method, several static and transient measurements were 
analyzed using physical models of the diode characteristics. To illus- 
trate the method, diffused silicon diodes were fabricated having substrate 

14 17 

(base) impurity concentrations ranging from 10 to nearly 10 phos- 
phorous atoms per cubic centimeter. The results show an emitter life- 
time that is much smaller than the base lifetime in the diode having 
the highest base doping concentration. In this diode, the recombination 
current from the emitter is 65% of the recombination current from the 
base, demonstrating the significance of the emitter in governing the 
static current-voltage dependence. The importance of emitter recombina- 
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tion to the transient characteristics was also demonstrated. Chap- 
ter III emphasized the techniques by which the base and emitter life- 
times are distinguished. It also demonstrated the need for carefully 
basing the quantitative analysis of the measurements on the underlying 
diode physics. 

The work of Chapter IV builds on the understanding and the pro- 
cedures developed in Chapter III. The methodology described in Chap- 
ter IV provides the first means available for the experimentally-based 
determination of the minority-carrier charge Q (effective gap shrinkage) 
and the effective lifetime in the emitter of p-n- junction solar cells. 

It also makes available for the first time a means for determining these 
material parameters in the emitter of any p-n- junction device. Its use 
was illustrated in the chapter by applying it to an n^-p solar cell of 
0.1 ohm-cm substrate resistivity. 

The methodology contrasts with the recent experimental study of 

bandgap shrinkage due to Slotboom and De Graaf [1] . That study restricted 

attention to gap shrinkage in the highly-doped base region of bipolar 

transistors. As a consequence the doping concentrations studied did 

not exceed 2x10 cn ; and p-type dopants only were considered in the 

particular study reported. In contrast the methodology described in 

Chapter IV determines the gap shrinkage in the emitter itself. It thus 

emphasizes the effects of the mechanical strain and of the very degene- 

19 -3 

rate doping concentrations » 10 cm ) occurring near the dif- 

fused emitter surface. These effects could be significant, particularly 
in shallow- junction emitters. In addition to the gap shrinkage, the 
methodology determines the effective carrier lifetimes in the base and 
the emitter. As is discussed in Chapter II, the emitter lifetime is 
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also likely to be sensitive to the effects of mechanical strain and 

degenerate doping concentrations occurring near the surface. 

Because of these differences, the methods of Slotboom and De- 

Graaff and those' presented here will 'likely give different results, 

and a comparison is in order. In their Eq. (18) , Slotboom and De Graaf 

give an empirical fit for gap shrinkage versus doping concentration. 

This fit holds for regions far' from the surface, for p-type dopants, 

19 -3 

and for doping concentrations less than 2x10 cm . But if the fit 
is nevertheless applied to the solar cell studied in Section IV, it 
yields the following values for the gap shrinkage in the shallow n-type 
emitter : 

20 —3 * 20 

AEg =125 meV for = 10 cm and AE^ = 135 meV- for = 2x10 cm 

These values are smaller than the gap shrinkage of 170 meV found from the 
methodology of Chapter IV. The approach of Slotboom and DeGraaff does 
not concern itself with lifetimes, which the methodology gives as 
Tg = 6.7visec and = 195 nsec for the particular cell studied. 

The methodology consists of a strategy for designing experiments 
and interpreting data consistently with the physical mechanisms governing 
device behavior. It is based on a model that has the following properties: 

(a) it describes both the static and the dynamic response; 

(b) it describes the device both in the dark and under 
illumination; 

(c) it characterizes the part of the device behavior coming 
from the surface and the bulk space-charge regions; 

(d) - it accounts for contributions to the device behavior 

coming from the -quasi-neutral emitter as well as from the 
quasi-neutral base. 
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As this model recognizes, many aspects of device behavior can come from 
processes not occurring within a single region of the device. But, 
for high-quality silicon solar cells,' the short-circuit current flowing 
in response to X-rays originates mainly in the quasi-neutral base region. 
Combining this response with the several other measurements indicated 
in Section 3.2 of Chapter IV fixes the value of the base lifetime x_ 
and- the related charge Q^q- The determination of .these parameters of 
the base helps toward the determination of lifetime and gap shrinkage in 
the quasi-neutral emitter. But the contribution to device behavior 
coming from the space-charge regions must be removed because it obscures 
that coming from the emitter. Several procedures are described for re- 
moving this contribution. They involve structural and electrical 
means, such as the use of an MOS guard ring, or heating the device, or 
resolving the data into components. The results of these procedures are 
then substituted into Eqs. (12) or (18) for the gap shrinkage in the 
emitter, and into Eqs. (13) or (19) for the emitter .lifetime. 

For the n -p cell structure studied in Chapter IV, the emitter 
plays a decisive role in determining solar-cell performance. The import- 
ance of the emitter has been suggested in earlier papers, but the exper- 
imental results given in Chapter IV provide direct evidence of this im- 
portance. For example, note that - 10. Thus the relaxation time 

B QN 

T measured by any of the methods discussed in Section 3.4 differs 
QN 

greatly from the effective base lifetime T^. . As was indicated in Chapter 
III, only for diodes with low base doping concentrations will 

B QN 

The treatment in Chapter IV has ignored the effects of areal in- 
homogeneity, which were described in Chapter II. These effects can be 
studied by applying the methods of this chapter to many small-area 
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device structures scribed from a large-area wafer subjected to solar- 
cell processing. 

. The methodology of Chapter IV bears on the' engineering design- 
of solar cells. It constitutes part of an effort to discover and 
describe the fundamental electronic mechanisms {Chapter II) influencing 
solar-cell efficiency, and to provide mathematical models that incorporate 
these mechanisms into solar-cell design. The methodology yields the 
excess minority- carrier charge stored in the emitter and the lifetimes 
in the emitter and the base. If the values of these parameters thus 
determined are correct, then any physical model that truly describes 
the cell must predict behavior consistent with these values. Because 
computer solution of the relevant differential equations is possible, 
physical models of any degree of complexity can be studied — models 
containing various combinations of the fundamental electronic mechanisms 
(Chapter II) that could influence device behavior. Comparison between 
the experimentally-based results of the methodology and the results 
predicted by computer simulation can then disclose which of the fun- 
damental mechanisms dominate in determining the cell efficiency. We 
have begun work aiming toward such a comparison. 

The work of Chapter IV relates mainly to the mechanisms of band- 
edge distortion and high recombination rates present in the diffused 
emitter. By providing methods for quantitatively determining these 
mechanisms, that work enables study of their relations to the processing 
used in solar-cell fabrication. It thus leads toward a systematic im- 
provement of open-circuit voltage and cell efficiency. A systematic 
study of the relations between processing and the electronic mechanisms 
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controlling efficiency is under way, and we plan to report the results 
of this study elsewhere. 
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APPLICATION OF THE SUPERPOSITION PRINCIPLE TO SOLAR-CELL 
ANALYSIS 
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List OjJ Symboti 

Electron and hole diffusivities 
Excess electron and hole densities 
Electric field 
Dielectric permittivity 

Optical generation rates of electrons and holes 
Solar cell current density, current 

Dark diode current density, current (subscript 0 denotes saturatic 

currei 

Electron and hole ciirrent densities 
Short-circuit current density, current 
Uncompensated photocurrent density 
Boltzmann ' s constant 

Electron and hole mobilities (b = u /y ) 

n p 

Electron and hole densities 

Ionized acceptor and donor impurity concentrations 
Intrinsic carrier density 
Negatively-charged trap density 
Electron charge 

Electron and hole steady-state recombination rates for 
thermal processes 

Series resistance 

Space charge density 

Absolute temperature 

Solar-cell terminal forward voltage 

Quasi-neutral base region voltage 

Intrinsic Fermi potential (or electrostatic potential) 

Junction space-charge region voltage 

Electron and hole quasi-Fermi potentials 
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I. INTRODUCTION 


Conventionally, the current in an illiiminated solar cell is de- 
scribed by the dark current- voltage characteristic Ij^ (V) shifted by 
the measured short-circuit current 

I(V) = -Id (V) t Ig^ (1) 

This description has served as a basis for many past theoretical esti- 
mates of various measures of solar cell performance. As is explained 
in Chapter I, the separation of I (V) into the components of (1) is use- 
ful because it essentially reduces the problem of explaining the 
discrepancy to the problem of understanding ^d ^^^ ’ dark current- 
voltage character is tics . 

As we show in this Appendix, however, the fundamental reasons 
justifying (1) hold only under restricted conditions. The use of (1) 
beyond its range of validity can lead to sizable errors in predicting 
the electrical characteristics of solar cells and thus to misleading 
conclusions about design. The shifting approximation holds, however, 
for the 0.1 ohm-cm silicon cells of interest in this research program. 

The appendix has several purposes : 

(a) To define the range of validity of (1) from funda- 
mental consideration of the relevant boundary-value 
problems ; 

(b) To note instances of practical importance for which 
the use of (1) will yield erroneous conclusions 
concerning cell performance; 

(c) To demonstrate, by laboratory experimentation and by 
exact computer solution of the pertinent boundary- value 
problems, the errors that can result from the use of (1) ; and 
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(d) To develop approximate analytic expressions for cell 
performance, based on fche device physics, that hold 
when (1) fails to apply. 

Despite the widespread use of (1) , it has apparently, until now, been 
justified only from an intuitive standpoint; no rigorous derivation 
of it has yet been done. In this Appendix, we provide such a derivation. 
The derivation is based on the principle of superposition , which applies 
if and only if the relevant boundary-value problems are linear. 

Following this derivation, we consider several examples of cell 
operation for which superposition and (l) do not apply. In these 
examples, we treat the effects of: (a) optical generation and net ther- 
mal recombination in the junction space-charge region; (b) high— injec- 
tion concentrations of holes and electrons; and (c) series resistance. 

The treatment of high injection involves consideration of cells in 
concentrated sunlight. Exact computer solutions of the boundary-value 
problems underlying cell operation together with laboratory ej^eriment- 
ation are used to support the theoretical developments. 


II. APPLICi\BILITY OF THE SUPERPOSITION 
PRINCIPLE TO SOLAR CELLS 

If a system is linear, its response to several excitations 
is the sum of the responses to each excitation applied alone. This is 
the principle of superposition , which is a basic theorem of linear dif- 
ferential equations [1] and which is widely used in such fields as 
electrical circuit theory [2] and control systems [3] . 
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we now appxy nne prxncipxe or superposirron cu a. JJ 711 — j mu- u-njn 
solar cell. As chapter II discusses, the fundamental equations describing 


this system in the dc steady state are [4,5]: 


J = +qD VN + qy NE 
N ^ n ^ n 


-qUnNVV^ 


j = -qD„VP + qy^PE 
P ^ P ^ P 


= -qypPWp 


0 = q(G^ - R^) 

0 = -VJp + q(Gp - Rp) 

V-E = -V^V^ =(p/e)= (q/e) “ N + - N^) 

J = + J . 

N P 


( 2 ) 

{ 2 a) 

(3) 
(3a) 

(4) 

(5) 

( 6 ) . 
(7) 


In general, this set of equations is nonlinear. Thus superposition 
does not apply unless approximations are made that remove the nonlin- 
earity from the system description. 

The first approximation to be made follows the widely-used Shock- 
ley method of analysis [4] which divides the cell into three coupled 
systems: the quasi-neutral emitter (QNE) and quasi-neutral base (QNB) 

regions, and the junction space-charge region (SCR) that separates them. 

The total current density J drawn from the cell results from the recom- 
bination-generation processes in the three regions. To demonstrate this 
for the one-dimensional n'^'-p solar cell of Fig. 1, we integrate the hole con- 
tinuity equation, (5), over the entire cell to obtain 


P B P E SCR 


^QNE 

+ q 


(Gp-Rp)dx +-q 


SCR 


(S-^) 


dx 


QNB 


,0 


(G„-Rp)dx . 


( 8 ) 


Because of the Ohmic contact present at the back surface (x = W ) , 


J = Jp(Wg) . 


(9) 
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The high recombination rates that occur in many cells near the 


illuminated surface (x = -W • - W ) imply that 

JS SCR 


J| » . 


P E SCR 
For the dc steady state, in general. 




^ 

Thus, using Eqs. (9), (10) and (11), we write (8) as 


J= q[ (G°-R)dx+ q 
'qne 


r 


SCR 


(G -R) dx+-q 

J 


(G°-R) dx 


QNB 


- ^QNE "^SCR "^gNB 


( 10 ) 


( 11 ) 


(12) 


Eq. (12) is analogous to Kirchoff s current law of electric 
circuit theory [2], as is illustrated in the equivalent-circuit re- 
presentation of Fig. 2. Eq. (12) and Fig, 2 are valid regard- ^ - 
less of whether or not superposition can be applied to the system 


describing the solar cell. 

We now examine the applicability of superposition to each of 
the regions of the cell- Consider first the quasi-neutral emitter 
and base regions. If low-injection conditions prevail in the quasi- 
neutral, regions and the effects of series resistance are negligible, 
then the expressions for the minority-carrier current and the net 


thermal recombination rate R both become linear in the minority- 
carrier concentration. Thus the continuity equation for minority 
carriers, (4), is linear, and the principle of superposition 
applies. The linearity of the system .holds in general despite the 
presence of built-in electric fields arising from impurity concen- 
tration gradients and despite heavy-doping effects [5] in the emitter 
such a’s energy-band-edge distortion. 
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To illustrate the use of superposition, consider the p-type 
base region of the solar cell of Fig. 1. For this system the inde- 
pendent excitations ar.e the optical generation rate G^(x) and the 
excess minority-carrier (electron) concentration AN(0). The res- 
ponse is the minority-carrier current Jgjjg (12) - Combining 
(11) and (12) , we write 


J = q 
QNB 


QNB 




(13) 


which is linearly related to the minority-carrier concentration 

and to the optical generation rate. Thus, using superposition, we obtain 


J = J 
QNB QNB 


+ J 




QNB 




|an(o) / 0 

The first term of (14) is a dark cxirrent -J 


AN(0) = 0 
D 


(14) 


QNB- 


The second 


term is a photocurrent J 


QNB- 


Therefore, 


,D 


■^QNB "^QNB "^ONB' 


(15) 


which can be represented by the equivalent circuit of Fig. 3. This 
decomposition of focuses attention on the dark current component, 

which is given by [4] 


QNB 


D , qVAT ' . 

= ^QBEO 


(16) 


for the assumed conditions of low injection and neglibible series 
resistance. 

Analogous results describe the current from the quasi- 

neutral emitter. Thus, if superposition applies, the component 
of current from the two quasi-neutral regions has the form of 
(1) ; it consists of a dark-current characteristic shifted by a 
photocurrent. 
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We now consider the applicability of superposition to the 
j'viuct-don space-chnrcis vociion, Tn this region the electti® 

E depends on the hole and electron concentrations, N and P. This ' 
introduces nonlinearity into the current density expressions, 

{ 2 ] and (3). Furthermore, the net thermal recombination rates, 
and , are nonlinear ly related to N and P. Therefore, the con- 
tinuity equations for holes and electrons, (4) and (5) , are 
both • nonlinear , and superposition does not apply to the junction 
space- charge region. 

The results of this section are illustrated by the equivalent- 
network representation of the solar cell in Fig. 4. The dashed 
lines that connect the current generators representing the quasi- 
neutral regions are closed or open depending on whether superpos- 
ition is applicable or not in the quasi-neutral regions. The breaks 
in connection in the space-charge region representation indicate that 
superposition is never valid in the junction space-charge region. 

Although the principle of superposition never applies to the 
junction space-charge region, the inapplicability does not necess- 
arily invalidate the resolution of the solar-cell current-voltage 
characteristics indicated in (1) . Using Fig. 4, we now consider 
special cases that illustrate this. 

2.1 Cases for which Eq. (1) applies 

Case A : Here the quasi-neutral regions dominate in determining 

both the photocurrent and the dark current. Low injection prevails 
throughout the quasi-neutral emitter and base, and series resistance ■ 
is negligible. Thus the dashed lines in Fig. 4 are short circuits-, 
indicating that superposition applies. From Fig. 4, 
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^ “ ^'^QNB "^gNE^ 

= “^gwR ~ - 1] 


(17) 


which is consistent with (i) _ 

Case B i Here the conditions are those stated for case A ex- 
cept that the space-charge region now also contributes significantly 
to the photocurrent. From Fig. 4, 


^ *'^gNR "^SCRj 




(18) 


Case C ; Here the dark current of the cell is dominated by the 
space-charge region and the photocurrent is dominated by the quasi- 
neutral regions. Series resistance is negligible. Thus, from Fig. 
4 and the theory of dark p-n junction diodes [8) / we have 


= ‘"gNR ■ 

Case D : Here the conditions are those stated in case A except 

that the space-charge and quasi-neutral regions both contribute 
significantly to the dark current. Thus, from Fig. 4 and ref- (6), 


J = J. 


QNR 




( 20 ) 


again in agreement with (1) . 

Discussion ; For cases A, B, and C, the open-circuit voltage 
has the simple form. 


V, 


OC 


mkT 

q 



( 21 ) 


where m == 1 for cases A and B, and m = m^^ for cases C. The open- 
circuit voltage for case D is the solution for V of (20) when 
J = 0, but no simple explicit expression for exists. 

The cases just described have practical importance. For ex- 


ample, cases A and B can describe a p-n-junction solar cell made 
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on a single-crystal silicon base having light or moderate doping 
concentration. Cases C and D can describe cells having highly- 
doped base regions. For such cells, the emitter region can contrib- 
ute significantly to the dark current [7] . 

To demonstrate the validity of (19) and (20) for cases C 
and D, we measured the illuminated and dark I-V characteristics of 

an n-on-p silicon solar cell having a highly-doped base region . 

17 -3 

(N. . = 10 cm ) . For this cell the thinness of the soace-charoe, 

= 

region causes the generation in the quasi-neutral base to dominate in determin 
the photocurrent. At voltages V near for one-sun illumination, the 
space-charge and quasi-neutral regions both contribute to the dark 
current (case D) . For the current range resulting from an illumination 
level of 1/100 s\in, -however, the recombination in the space-charge 
region alone determines the dark current (case C) . Our measured data 
showed excellent agreement with (1) for both levels of illumination. 

Further experimental evidence of the validity of (1) for the low- 

■ resistivity silicon cells of interest in this report appears in Chapter 

IV where we found that the value for the current determined by 

QNO 

^^J^k I-V measurements agreed well with that determined by illuminated 
I-V measurements. 

2.2 Cases for which Eg. (1) does not apply 

Any of the following conditions will invalidate the principle of 
superposition and thus (1): 

(a) the space-charge region contributes importantly to both 
the photocurrent and the dark current; or 

(b) the carrier concentrations in the quasi-neutral regions 
exceed low- injection levels; or 
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(c) the series resistance contributes significantly to the cell 
current- volt age characteristics; or 

(d) the material parameters, such as the minority- carrier life- 
time, depend on the illumination level; or 

(e) the volume of the region producing the photocurrent changes 
appreciably as the cell is loaded (and the terminal voltage 
V changes) . 

The invalidity of (1) for CdS cells has been discussed and shown 
experimentally [8] . The invalidity of (1) for amorphous silicon (a-Si) 
cells [9] and the consequent poor fill factor probably arises, in part, 
from the decrease in vol\ame of the junction space-charge region that 
accompanies loading of the cell. Because the diffusivity and hence 
the diffusion length of the quasi-neutral regions of a-Si cells is 
small, the photocarriers generated in the junction space-charge region 
dominate in determining the photocurrent I. Thus the invalidity of (1) 
for a-Si cells owes its origin in part to condition (e) listed above and, 
in part, to condition (a). 

Experiments showing the invalidity of (1) for single-crystal silicon 
cells will be discussed in the following sections, which treat the solar- 
cell characteristics for cases (b) and (c) . 
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III. HIGH INJECTION IN THE QUASI-NEUTRAL BASE 


Ah important example of high injection occurs in the guasi-neutral 
base region of a solar cell exposed to concentrated sunlight [10] . Despite 
the inapplicability of superposition to this cell from a strict stand- 
point,- an analytic description- based on superposition can be developed, 
for the case of extreme high-level injection where P==N in the base re- 
gion. For this case, superposition yields [103 


J = J (Q) = fO) 
N N 


0 

AN(0) ^ 0 


t J^(0) 


G° 5^ 0 
AN(0) = 0 


= ^ -^UPC • ^22)- 

here. In (22) , is the dark current-voltage characteristic of- 

the cell. The teacn in (22) is the photocurrent in a highly 

illuminated cell to which a reverse bias is applied 

to insure AN (0) = 0; is essentially»the maximum current that can 

be drawn from the solar cell for a given level of illumination. We thus 


call it the uncompensated photocurrent . In general, J 

Using (22) and including effects of high injection in the 
base region, one can derive analytical expressions for the open-circuit 
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voltage of highly-illuminated silicon solar cells [10] . In Fig. 5, 

-f 

values of calculated from these expressions for lOfi-cm n -p cells 
[(V proportional to (2/b+l) (kT/q)ln{J. /J )] and 4fi-cm p^-n cells 
'[{V^^ proportional to (2b/b+l) (kT/q) ln(J^^/j^^) ] are plotted against 

These dependencies of on compared at several levels of 

illuminatxon with values of V calculated from exact nimerical solutions 

OC 

of (2) through (7) - These results are contrasted with the classical 
relation, derived from (1) , 

- (kT/gjlnW^c/V ' 

which is a special case of (21) . At low-injection levels, J = ^cn' 
otherwise, J a considerable potential difference can 

develop across the quasi-neutral base [10] . The inapplicability of 
(1) is apparent in Fig. 5. This results because the system is not 
strictly linear [10] . 

The theoretical results are supported by the experimental data 
plotted in normalized fashion against the measured three tj^es 

of n^-p solar cells: a cell having a base resistivity of 10 ohm-cm, one . 

made with a back-surface field (BSF) and the other without a BSF. At 
high levels of illumination (large / high injection occurs in the 

base of the 10 ohm-cm cells at the open-circuit condition. For the 10- 
ohm-cra cell without a BSF, the measured dependence of on J^p^^ agrees 

with 'the theoretical' predictions of Fig. 5 and demonstrates experimentally 
that (1) is invalid. The dependence observed for the BSF cell again 
demonstrates the inapplicability of superposition. The nonlinearities 
associated with high injection in the base region of this cell present 
•an analytical problem different from that describing a non-BSF cell. 
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Basic corrections of classical theory required by these npnlinearities 

will be treated in a later publication. 

High injection does not occur in the 0.3fJ-cm cell, and its 

V vs. J characteristic agrees with the classical result expressed 
OC UPC 

in (23) - 

IV. SERIES RESISTANCE . 

In a p-n- junction solar cell, series resistance R can originate 
in the metal contacts, in the metal- semiconductor interface, in the 
emitter, and in the base. The presence of R invalidates the direct use 

O 

of superposition. To illustrate this, consider, as the system under 
study, a p-type quasi-neutral base region. Though the electron contin- 
uity equation, (4) , remains linear for low injection, the presence 
of Rg results in a nonlinear boundary condition relating An and J^: 

AN(0) “ [exp(qV/kT)exp(qIR /kT) - 1] , (24) 

where I depends on J (0) . The nonlinearity of (24) precludes the 

N 

use of superposition in which AN(0), defined by (24), is regarded 
as an independent excitation. 

A solar cell with significant R can be treated analytically by 
separating R„ from the intrinsic system, as in Fig. 7. The intrinsic 
system is defined so that 

AN(0) [exp(qV /kT) - 1] . (25) 

IS 

Since AN(0) thus defined depends only on the voltage V^g across the systera, 
it can be regarded as an independent excitation; Thus the intrinsic 
system of Fig. 7 can be treated as described in Section II, where we 
indicated the conditions for which superposition applies - 
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Use of superposition requires that An'(O) be set to zerO/ which de- 
fines the photocurrent As for the high- inject ion condition des- 

cribed in Section III., here also a short circuit applied across the ceil 
terminals -fails to make AN (0) because 


V = V 

IS 


I R_ 


(26) 


From {25) we see that AN(0) = 0 implies V = 0, and from (26) 

-i-O 

that a reverse bias of about I R must therefore be applied to the - 
solar cell to experimentally determine We note again that J^p^^ 

> Jgj., where = J (v=0) . 

For low- injection conditions, R is a constant, independent of I, 
and can be determined eiqierimentally [H]. 

V. DISCUSSION 

The- basis underlying this Appendix is the superposition principle. 

Its use here provides insight about the approximation of (1) that 
the characteristic of an ill\nninated solar cell is the dark current- 
voltage characteristic shifted by the short-circuit current. In- 
terpreted from the standpoint of the superposition principle, the current 
I(V) in ( 1 ) is the sum of the response to two excitations applied 
independently, one at a time. One excitation is the hole-electron op-, 
tical generation rate resulting from the sunlight. The other is the 
excess carrier concentration at each edge of the junction space-charge 
region, which is related to the voltage difference across the cell .terminals* 
The two responses are therefore the short-circuit current and the dark- 
current characteristic. ' 
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The optical generation rate and the excess carrier concentration 
at the SCR edges are chosen as the excitations because they are the 
independent variables in the boundary-value problem describing the cell. 

These variables are linearly related to the current and to the excess carrier 
concentrations within the cell under certain well-defined conditions 
which, therefore, are necessary for superposition to apply. These 
conditions are : 

(a) the junction space-charge region may contribute importantly 
to either the photocurrent or the dark current, but not 

to both; 

(b) the -carrier concentrations in the quasi-neutral regions 
must stay within low-injection levels; 

(c) the series resistance must contribute negligibly to the 
cell current— voltage characteristics; 

(d) the material parameters, such as the minority-carrier 
■lifetime, must be essentially independent of the illumi- 
nation level . and 

(e) the volume of the regions that contribute appreciably 
to the photocurrent must stay essentially constant as 
the cell is loaded. 

We have described corrections to superposition-based 
theory if conditions (b) or (c) are violated. The violation of condi- 
tion (b) discussed here relates to an application of much practical importance — 
the performance of (single-crystal) silicon p-n junction cells in concen- 
trated sunlight. Although in a strict sense the superposition principle 
is invalid for this application, we used superposition in developing the 
correction. This use led to a theory that agreed well both with labor- 
atory experimentation and with the exact computer solution of the underlying 
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differential equations (2) - (7) . 

This is the first application of the superposition principle of 

differential-equation -theory to the analysis of solar cells. We be- 
lieve that this principle offers valuable insight about the range of 
validity of (1) and of the considerable simplifications that (1) yields. 
Beyond -this, the treatment discussed here of a cell in concentrated 
sunlight has shown that superposition provides a basis from which to 
develop corrections when (1) fails to hold. We anticipate that the 
superposition principle will prove useful also in applications other 
than those discussed in detail here, including applications to cells 
made of materials other than single-cirystal silicon. 


113 



REFERENCES FOR APPENDIX 


1. P. w. Berg and J. L. McGregor , Elementary Partial Differential Equa- 
tions , p. 7, Holden-Day, San Francisco, 1963. 

G. A. Korn and T. M. Korn, Mathematical Handbook for Scientists and 
Engineers , Sections 9.3-1, 10.4-1, 10.4-2, and 15.4-2, McGraw-Hill, 1961. 

2. D-F. Tuttle, Network Synthesis , p. 19, John Wiley & Sons, 1958. 

II. H. Skilling, Electrical Engineering Circuits , p. 353, ®John Wiley 
S Sons, 1965. 

3. J. L. Bower and P. M. Schultheiss, Introduction -do the Design of 
Servomechanisms , pp. 11, 100, John Wiley S Sons, 1958. 

4. w. Shockley, "The Theory of p-n Junctions in Semiconductors and p-n 
Junction Transistors," Bell Sys, Tech. J. , Vol. 28, pp. 435-489, 

July 1949 . 

•5., F. A. Lindholm and.C. T. Sah, "Fundamental Electronic Mechanisms 
Limiting the Performance of Solar Cells," IEEE Trans. Electron 
Devices , April 1977. 

6. C. T. Sah, R. N. Noyce, and W. Shockley, "Carrier Generation and 
Recombination in p-n Junctions and p-n Junction Characteristics," 

Proc. IRE , Vol. 45, pp. 1228-1243, Sept. 1957. 

C. T. Sah, "Effect of Surface Recombination and Channel on p-n 
Junction and Transistor Characteristics," IRE Trans, on Electron 
Devices , Vol. ED-9, pp. 94-108, Jan. 1962. 

7. F. A.. Lindholm, A, Neugroschel, C.T. Sah, M. P. Godlewski, and H. W. 
Brandhorst, Jr., "Methodology for the Experimental Detemination- of 
Gap Shrinkage and Lifetimes in the Emitter and Base of Solar Cells and 
other p-n Junction Devices," IEEE Trans. Electron Devices , April 1977; 
also Record of Twelfth Photovoltaic Specialists Conf . , Nov. 1976. 

•8. 'K-. W. Boer, "CdS/Cu S Hetero junction Model," Proc. Sharing the Siin , 

Int. Conf., Photovoltaics , Vol. 7, pp. 130-145, Winnipeg, Manitoba, 
Canada, Aug. 1976. 

9. C. R. Wronski, D. E. Carlson, R. E. Daniel and R. R. Triano, "Electri- 
cal Properties of a-Si Solar Cells," Technical Digest 1976 Int. Electron 
Devices Meeting (76CH1151-OED) , pp. 75-78, Dec. 1976; also D.E. Carlson 
and C.R. Wronski, Appl. Phys. Lett ., vol. 28, No. 11, P. 602 (1976). 

10. J. G. Fossum and F. A. Iiindholra, "The Dependence of Open-Circuit Voltage 
on Illumination Level in p-n Junction Solar Cells," IEEE Trans. Electron 
Devices , April 1977 . 

11. M. Wolf and H. Rauschenbach, "Series Resistance Effects on Solar Cell 
Measurements," Advanced Energy Conversion , Vol. 3, pp. 455-479, Apr.- 
June, 1963. 


114 




»gS«VW?iVi^. 






Ids regardless of the applic- 


L6 




Fig. 3 Equivalent circuit for the quasi-neutral base if super- 
position applies. A siniilar equivalent circuit represents 
the quasi-neutral emitter - 
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(a) Entire solar cell in which AN{0) depends nonlinearly on I; 

(b) Intrinsic solar cell in which AN (0) is proportional to 
[exp(q^^^T) -1], 
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